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ABSTRACT 

A technique for the use of a digital computer for simul¬ 
taneous generation of feedback control system transient 
response and sensitivity to parameter variations is developed. 
The sensitivity model for a basic linear control system is a 
cascaded replica of the original system with the sensitivity 
pickoff points corresponding to the error and feedback signals 
of the original system. This is convenient for computer anal¬ 
ysis with DSL/360. For systems with nonlinear elements, the 
sensitivity model is implemented by making the instantaneous 
gain equal to the slope of the nonlinear function. Sensi¬ 
tivity analysis is utilized for improvement of system 
performance, determining parameter tolerances, and predicting 
changes in the response due to parameter variations. 
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INTRODUCTION 


I. 


The early development of sensitivity analysis was 
applied to circuit theory where it is useful to know how 
sensitive the response of the circuit is to variations of 
circuit parameters such as resistance, capacitance,, and 
inductance. There are many different definitions of sensi¬ 
tivity available. H. Bode [1] defines a logarithmic or 
normalized sensitivity function of a variable, x^(t,a^), 


a £n a . 


and I. M. 



3£n 


a 


j 


Hortowitz [2] adapted the inverse of this definition which 
is widely used. J. Leeds [3] presented sensitivity function 
for incremental changes in parameters and S. R. Parker [4] 
recently generalized the technique of sensitivity analysis for 
nonlinear circuits. Much of this work entailed the definition 

of R. Tomovic [5] who defined the sensitivity coefficients as 

X. 3x. X. 3x. 

1 _ i .11 

= K— . Kokotovic [6] adapted the definition S = tto- 

aj 3aj i J f 

for use in control systems and developed a method of generat¬ 
ing sensitivity coefficients using analog computers. In this 
thesis the method of Kokotovic [61 is expanded for use with 
digital computer techniques utilizing block diagxams and 
transfer functions of forward and feedback loops in control 
systems. The IBM/360 Digital Simulation Language (DSL) 
program [7] is very well suited for modeling a control system 
and obtaining the time response and sensitivity coefficients 
simultaneously. 
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DSL is a continuous system simulator whose non-procedural, 
problem-oriented input language has the attractive building- 
block approach of digital-analog simulators with the added 
power of logical and algebraic equation notation. The DSL 
system provides a basic set of 35 blocks from which a phys¬ 
ical system may be modeled. Included in this set are such 
conventional analog components as integrators, comparators, 
and transfer functions. The Kokotovic [6] method of sensi¬ 
tivity analysis is readily adaptable to the DSL program as 
control systems in block-diagram form can be programed in 
the s-domain. 

Section II of this thesis applies the Kokotovic [6] 
method of generating sensitivity coefficients of the basic 
feedback control system for the forward and feedback gain 
parameters. Sensitivity functions are derived and applied 
to a third-order system. Examples are given for the use 
of sensitivity coefficients in analysis, parameter adjust¬ 
ment for improvement of system performance, and parameter 
tolerance considerations. Section III applies the 
Kokotovic [6] method of sensitivity analysis to multiloop 
systems. Gain and time-constant sensitivity coefficients 
are derived and applied to a regulated power generator as an 
example. Section IV demonstrates the Kokotovic [6] method 
for the general case of differential equations, and Section 
V is devoted to the nonlinear system. 

Sensitivity analysis of a control system offers an 
analytical method for determining those parameters of a 
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dynamic system which have the greatest influence on the 
system performance. High-speed digital computers give 
very precise solutions for the system equations, but it is 
impossible to realize exact mathematical models and para¬ 
meter values. Additional variations in parameters may also 
be caused by wear, aging, shifting of operating point, 
temperature, and etc. It is beneficial to know the effect 
of these small parameter variations on system performance. 
Sensitivity analysis can give this necessary insight so 
that all parameters of a dynamic system are not given equal 
importance. This is more efficient than repetitive computer 
solutions with parameters changed by increments. Analysis 
utilizing sensitivity coefficients can put emphasis on 
specific component tolerances and points that are critical 
to the performance specifications. 
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II. BASIC CONTROL SYSTEM SENSITIVITY 


A. INTRODUCTION 

The Kokotovic [6] method for computing sensitivity 
coefficients for various parameters is restricted to linear, 
time-invariant, deterministic systems. Although it is 
true that physical systems do not have these ideal properties, 
linear approximations for these systems have proven adequate 
in many useful cases. When parameter values are to be 
approximated, it is important to estimate the influence of 
each parameter on the system response. This can be done 
accurately by sensitivity analysis. Also, limits of indi¬ 
vidual parameter variations may be deduced and tolerance 
levels formulated. 

The Kokotovic [6] method of solving for sensitivity 
coefficients has the advantage that all of the functions 
are generated simultaneously with the time-domain solution. 
Kokotovic [6] gives analytic proof that this is possible. 

As an engineering approach to analysis, it is common to 
approximate a high-order system by an equivalent second- 
order system to enable using graphs and tables for second- 
order response characteristics. Qualitative observations 
of graphs of sensitivity coefficients can give insight into 
the transient performance in terms of overshoot, settling 
time, rise time, and etc., for higher degree characteristic 
equations as will be seen in examples. 
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B. GAIN CONSTANT SENSITIVITY COEFFICIENTS 


The sensitivity coefficient S„ is defined as the ratio 
of the change in output to the fractional change in the 
parameter K. 


_ 3C _ 3C 
K “ 3K/K 3 5,nK 


( 2 . 1 ) 


where C is the controlled output and K is the parameter being 
investigated. 

The canonical form of a feedback control system after 
Laplace transformation is shown in Fig. 2.1. The definitions 
of terms are as follows: 


K 

KG(s) 

F 

FH (s) 
KG(s)FH(s) 
R(s) 
E 
B 


direct gain constant 
direct transfer function 
feedback gain constant 
feedback transfer function 
loop transfer function 
reference input 
error signal 
feedback signal 


The basic equations are as follows: 


^ KG(s) _ 

R ^ l+KG(s)FH(s) 


_ 1 

R^®^ ~ l+KG(s)FH(s) 


B. ^ KG(s)FH (s) 

R^ ^ l+KG(s)FH(s) 


(2.2a) 

(2.2b) 


(2.2c) 


In the work that follows, capital letters denote the 
Laplace transformation of the functions described except for 
the letters K and F which are constants. 
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R^ 



KG (s) 

C 

- 

1 






B 

FH(s) 






Fig.2.1. Canonical form of a feedback, control system. 



Fig.2.2. Feedback control system and sensitivity model. 
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TIti£(SEC) 


Fi6.2.3. 


Basic control systero response and sensitivity 
for K = 50 . 
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Applying (2.1) to (2.2a) yields the gain sensitivity 
coefficients 


RKG 


C 


S 


K 


(1+KGFH)^ 


1+KGFH 


(2.3) 


S 


-RK^G^FH 


C(-KGFH) 

1+KGFH 


(2.4) 


F 


(1+KGFH)^ 


These two gain sensitivity functions are shown in block 
diagram in Fig. 2.2. This shows that the sensitivity model 
is an exact replica of the original system with the controlled 
output as the input to the sensitivity model. The error and 
feedback signals' of the sensitivity model are the desired 
sensitivity functions. It should be noted that the sensi¬ 
tivity model requires no new components or transfer functions 
and all quantities are generated simultaneously. 

C. APPLICATION 

In this thesis the equations for the feedback systems 
are programed on the IBM 360 computer using the Digital 
Simulation Language (DSL) program for both the transient 
solution and the sensitivity coefficients. This program 
has the advantages of programing in the s-domain in trans¬ 
fer function form and the accuracy of the digital computer 
with graphical and numerical results. A sample program is 
included in Appendix A and the numerical solution and 
sensitivity coefficients are listed in Table 2.1. Fig. 2.3 
is the graph of the output and sensitivity coefficients. 
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The direct transfer function is 


KG 


50 ^ 

s (s+4) (s+6) 


(2.5) 


For example/ data from Table 2.1 when is a maximum 

Time T = 1.0 

Gain K = 50 

Sensitivity S = 0.70863 

K 

Input R = unit-step 

Output C = 0.96456 

For large variations (2.1) becomes 

SK = ^ (2.6) 

Let the gain, K, change 10%, then 

ak 

AC - Sr ^ = 0.70863(0.1) = 0.070863 


This shows that the output will change approximately + 0.071 

if the gain varies from 50 to 55 or 45. This can be 

verified from Table 2.2 for the following data: 

K = 55 

T = 1.0 

C = 1.0331 

C = C + AC = 0.96456 + 0.070863 = 1.03541 

and for 

K = 45 

T = 1.0 

C = 0.89133 

C = C-AC = 0.96456 - 0.07086 = 0.89370 
This shows that sensitivity coefficients for changes of gain 
as large as 10% produce accuracy within 0.3% when compared to 
incremental intervals as denoted in partial derivatives. 
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Another example for the use of sensitivity coefficients 
is to compute the gain required to reduce the peak overshoot 
of the system transient response to a more desirable specifi¬ 
cation of one-half its initial value. Table 2.1 gives the 
following data: 


K = 50 

Peak overshoot = 0.2 

T = 1.5 

= 0.36568 

Rearranging (2.6) 


AK 


ACK 

Sr 


K' 


K - AK 
50 (1 


K(l—) 

0.1 

0.36568 


) 


36 


(2.7) 

( 2 . 8 ) 


Table 2.2 for K = 36 shows that the peak overshoot is reduced 
to 0.1 and Fig. 2.4 shows the reduction of peak overshoot by 
a factor of one-half. 

Any other design specification, such as settling time or 
rise time, can be adjusted by utilization of the sensitivity 
coefficient. In the example it should be noted that reducing 
the overshoot by reduction in gain, increases the rise time 
by 0.4 second. 

Tolerance considerations for parameters are encompassed 
within the definition of sensitivity coefficient. 


Rearranging (2.6) 

Ar 

Tolerance = AK/K = — 


(2.9) 
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1.2 


1.0 


0.8 


0.6 


C.4 


0.2 


C.O 



O.C 2.0 4.0 6.0 8.0 10.0 

TIME(SEC) 


Fig.2,4. Overshoot reduction of 0.1 hy utilization of 
sensitivity coefficients. 
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0.0 2.0 4.0 6,0 8.0 10. 

TIME (SEC ) 


Fig.2.5. Gain tolerance of l4^ for overshoot of i0.05. 
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BASIC SERVO SENSITIVITY COEFFICIENTS 


TIME 

OUT 

SENSK 

0.0 


0.0 

0.0 

l.OOOOE- 

■01 

6.5182E-03 

6.516CE-C3 

2.00C0E- 

•01 

4.1385E-02 

4.1258E-02 

3.0000E- 

■01 

1.1165E-01 

1.1054E-ri 

4.0000E- 

■01 

2.1292E-01 

2.0812E-C1 

5.0000E- 

■01 

3.3623E-01 

3.2212E-C1 

6.0000E- 

■01 

4.7147E-01 

4.3885E-C1 

7.0000E- 

■01 

6.0911E-01 

5.4526F-C1 

8.00C0E- 

■01 

7.41C5E-C1 

6.3042E-01 

9.0000 E- 

■01 

8.6100E-01 

6.864CE-r 1 

l.OOOOE 

00 

9.6456E-01 

7.0863E-C1 

I.IOOOE 

00 

1.0491E 00 

6.959CE-C1 

1.2000E 

00 

1.1137E 00 

6.5CC5E-01 

1.3000E 

00 

1.1586E 00 

5.7542E-C1 

1.40C0E 

00 

1.1853E OC 

4.7816E-01 

1.5000E 

00 

1.1959E 00 

3.6558E-01 

1.6000E 

00 

1.1931E 00 

2.4539E-01 

1.7000E 

00 

1.1797E 00 

1.25 If E-01 

1.8000E 

00 

1.1587E 00 

1. 1496E-C2 

1 .9000E 

00 

1.1328E 00 

-8.9767E-02 

2.0000E 

00 

1.1046E 00 

-1.7441E-01 

2.lOOOE 

00 

1.076CE 00 

-2.3964E-01 

2.2000E 

00 

1.0489E 00 

-2.8418F-C1 

2.3nC0E 

00 

1.0245E 00 

-3.0811E-ri 

2.A000E 

00 

1.0036E 00 

-3.1273E-C1 

2.5000E 

00 

9.8672E-01 

-3.0032E-01 

2.6000E 

CO 

9.7404E-01 

-2.7383E-01 

2.7000E 

GO 

9.6544E-01 

-2.3663E-C1 

2.8O00E 

00 

9.6056E-01 

-1.9225E-C1 

2.900CE 

00 

9.5895E-01 

-1.4414E-C1 

3.0000E 

00 

9.6002E-01 

-9.5469E-f2 

3.irCOE 

00 

9.6318E-01 

-4.8996F-C2 

3.2000E 

00 

9.6782E-01 

-6.9429E-03 

3.300CE 

00 

9.7336E-01 

2.9C47E-02 

3.4000E 

00 

9.7932E-01 

5.7914E-C2 

3.5C00E 

00 

9.8525E-01 

7.9l65E-f2 

3.60CCE 

00 

9.9r83E-01 

9.2808E-C2 

3.70n0E 

00 

9.9581E-01 

9.9281E-C2 

3.8000E 

00 

i.oonoE 00 

9.9353F-C2 

3.9C00E 

00 

1.0034E 00 

9.4037E-C2 

4.0000E 

00 

1.0C58E 00 

8.4A86F-C2 

4.1000E 

00 

1.0075E 00 

7.1906F-C2 

4.2000E 

00 

1.0084E 00 

5. 7483E-02 

4.3000E 

00 

1.0086E 00 

4.2310E-C2 

4.4000E 

00 

1.0083E 00 

2.735CE-C2 

4.5000E 

00 

1.0075E 00 

1.3390 E-C2 

4.6000E 

00 

1.0065E 00 

1.0386E-03 

4.7000E 

00 

1.0053E CO 

-9.2916E-03 

4.8000E 

CO 

1.0041E 00 

-1.7368E-C2 

4.9000E 

00 

1.0n28E 00 

-2.3124E-r2 

5.0000E 

00 

1.0017E 00 

-2.6629E-C2 

s.ioncE 

00 

1.0007E 00 

-2.8C75F-02 

5.2000E 

00 

9.9984E-01 

-2.7730E-C2 

5.3000E 

00 

9.9918E-01 

-2.5919F-C2 

5.4000E 

00 

9.9869E-01 

-2.2996E-02 

Table 2.1. Computer output data for unity 




SENSE 

- 0.0 

-2.2129E-06 
-1.2685E-0A 
-1.1099E-03 
-4.7931E-03 
-1.4109E-C2 
-3.2619E-C2 
-6.3846E-C2 
-1.1063E-01 
-1.7460E-01 
-2.5593E-ni 
-3. 5322E-01 
-4.6363E-01 
-5.8319E-01 
-7.0712E-C1 
-8.3030E-C1 
-9.4767E-01 
-1.0546E OC 
-1.1472F 00 
-1.2226E 00 
-1.2790E or 
-1.3157E 00 
-1.3331E OC 
-1.3326E or 
-1.3163E on 
-1.2870E OC 
-1.2479E 00 
-1.2021E 00 
-1. 1528E 00 
-1.1031E or 
-1.0555E OC 
-1.0122E 00 
-9.7476E-01 
-9.4432E-01 
-9.2140E-C1 
-9.0609E-01 
-8.9802E-01 
-8.9653E-01 
-9.0066E-01 
-9.0933E-C1 
-9.2136E-01 
-9.3558E-01 
-9.5088E-01 
-9. 6627E-C1 
-9.8091E-01 
-9,9414E-C1 
-1.0055E 00 
-1.0146E OC 
-1.0214E 00 
-1.0260E OC 
-1.0283E or 
-1.0288E OO 
-1.0276E no 
-1.0251E 
-1.0217E or 


feedback system. 
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BASIC SERVO SENSITIVITY COEFFICIENTS 


Table 2 


TIME 

s.snooE 00 
5.6000E 00 
5.7COOE 00 
5.B000E 00 
5.9000E 00 
6.0000E 00 
6.1000E 00 
6.2000E 00 
6.3000E 00 
6.4000E 00 
6.5C00E 00 
6.6000E 00 
6.7000E CO 
6.8000E 00 
6.9C00E 00 
7.0000E 00 
7.1000E 00 
7.2000E 00 
7.3000E 00 
7.4000E no 
7.500CE 00 
7.6000E 00 
7.7COOE 00 
7.8000E 00 
7.900CE 00 
8.0000E 00 
8.100CE 00 
8.2000E 00 
8.3C00E 00 
8.40C0E 00 
8.5000E 00 
8.60CCE 00 
8.70Q0E 00 
8.8000E 00 
8.9000E 00 
9.0000E 00 
9.1000E 00 
9.2000E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000E 00 
9.7C00E 00 
9,8n00E 00 
9.900CE 00 
l.OOOOE 01 


OUT 

9.9838E-01 
9.9823E-01 
9.9821E-01 
9,9829E-0l 
9.9846E-01 
9.9868E-01 
9.9893E-01 
9.9920E-01 
9.9945E-01 
9.9969E-01 
9.9989E-01 
I.OOOIE 00 
1.0002E 00 
1.0003E 00 
1.0003E 00 
1.0C04E 00 
1.0004E 00 
1.0003E on 
1.0003E 00 
1.0003E 00 
1.0002E 00 
1.0002E 00 
I.OOOIE 00 
I.OOOIE 00 
l.OOOOE 00 
9.9998E-01 
9.9995E-01 
9.9994E-01 
9.9992E-01 
9.9992E-01 
9.9992E-01 
9.9993E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9997E-01 
9.9998E-01 
9.9999E-31 
9.9999E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 


-1.9314E-02 
-1.52C8E-02 
-1.0980E-02 
-6.8853E-C3 
-3.1281E-r3 
1.4061E-04 
2.8257E-03 
4.8800E-03 
6.3002E-C3 
7.1194E-C3 
7.3986E-C3 
7.2171E-03 
6.6686E-03 
5.8481E-03 
4.8504E-03 
3.7631F-03 
2.6627E-C3 
1.6125E-f3 
6.6310E-04 
-1.5068E-r4 
-8.0872E-04 
-1.3018E-C3 
-1.6327E-C3 
-1.8120E-f3 
-1.8587E-03 
-1.7906E-C3 
-1.6360E-03 
-1.4186E-C3 
-1.1616E-C3 
-8.8680E-C4 
-6.1244E-C4 
-3.5483E-C4 
-1.2463E-C4 
6.9797E-05 
2.2489E-C4 
3.3873E-C4 
4.1258E-C4 
4.4990E-C4 
4.5526E-04 
4.34C4E-04 
3.9238E-r4 
3.3724E-04 
2.7323E-04 
2.0635F-04 
1.4C67E-04 
8.0466F-C5 


SENSF 

-1.0177E 00 
-1.0134E 00 
-1.0092E 00 
-1.0052E 00 
-1.0016E or 
-9.9854E-C1 
-9.9611E-01 
-9.9432E-01 
-9.9315F-01 
-9.9257E-01 
-9.9249E-01 
-9. 9284E-01 
-9.9352E-01 
-9.9444E-01 
-9.9549E-01 
-9.9661E-01 
-9.9771E-01 
-9.9873E-C1 
-9.9965E-01 
-1.0004E 00 
-l.OOlOE 00 
-1.0015E 00 
-1.0017E 00 
-1.0019E CC 
-1.0019E OC 
-1.0C18E CC 
-1.0016E Of 
-1.0014E 00 
-l.OOllE 00 
-1.0008E 00 
-1.0005E 00 
-1.0C03E OC 
-I.OOOIE 00 
-9.9987E-01 
-9.9973E-C1 
-9.9963E-01 
-9.9956E-01 
-9.9954E-C1 
-9.9954E-C1 
-9.9957E-01 
-9.9961E-01 
-9.9967E-01 
-9.9974E-01 
-9.9981E-01 
-9.9987E-01 
-9.9993E-01 


.1. (Continued) 
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TIMF 

0.0 

l.OOOOE-Cl 

2.0000E-01 

3.0000E-01 

4,0000E-ni 

5.0000E-01 

6.0000E-01 

7.0CC0E-01 

e.ooooE-oi 

9.0CC0E-01 
l.OOOOE CO 
I.ICOOE 00 
1.2C00E 00 
1.3000E 00 

1. AOOOE 00 
1,5000E 00 
1.60C0E 00 
1.7000E 00 
l.SOOCE 00 
1.9C00E 00 
2.0000E 00 
2.1000E 00 
2.2000E 00 
2.3000E 00 

2. A0C0E 00 
2.50C0E 00 
2.6COOE 00 
2.7000E 00 
2.8000E 00 
2.9000E 00 
3.0000E 00 
3.1000E 00 
3.2000E 00 
3.3000E 00 

3. A0C0E 00 
3.5000E 00 
3.6000E 00 
3.7CC0E 00 
3.8r00E 00 
3.9000E OC 
A.OOOOE 00 
A.IOOOE 00 
A.2000E 00 
A.3000E 00 
A.AOOOE 00 
A.5000E 00 
A.6000E 00 
A.7000E 00 
A,8000E 00 
A.9000E 00 
5.0000E 00 
5.1C00E 00 
5.2000E 00 
5.300CE 00 
5.A000E 00 


Table 2.2. 


K =50 K=55 K;=45 


K=36 


C.O 

6.5182E-03 
A.1385E-02 
1.1165E-01 
2.1292E-01 
3.3623E-01 
A.71A7E-ni 
6.091lE-01 
7.A105E-01 
8.6100E-01 
9.6A56E-01 
1.0A91E 30 
1.1137E 00 
1.1586E 00 
1.1853E 00 
1.1959E 00 
1.1931E 00 
1.1797E 00 
1.1587E 00 
1.1328E 00 
1.10A6E 00 
1.0760E 00 
1.0A89E 00 
1.02A5E 00 
1.0036E 00 
9.8672E-01 
9.7A04E-01 
9.65AAE-01 
9.6056E-01 
9.5895E-01 
9.6002E-01 
9.6318E-01 
9.6782E-ni 
9.7336E-01 
9.7932E-01 
9.8525E-01 
9.9083E-01 
9.9581E-01 
l.OOOOE 00 
1.003AE 00 
1.0058E 00 
1.0C75E 00 
1.008AE 00 
1.0C86E 00 
1.0083E 00 
1.0075E 00 
1.0065E 00 
1.0053E 00 
I.OOAIE 00 
1.0028E 00 
1.0017E 00 
1.0007E 30 
9.998AE-01 
9.9918E-01 
9.9869E-01 


OUT 

0.0 

7.1698E-03 
A.5510E-02 
1.2269E-31 
2.3363E-01 
3.6830E-01 
5.1503E-01 
6.63C1E-01 
8.0304E-01 
9.2801E-01 
1.0331E 00 
1.1156E 00 
1.17A7E OO 
1.211AE 00 
1.2276E 00 
1.226AE 00 
1.2112E 30 
1.1857E 00 
1.1536E 00 
1.1180E 00 
1.0820E OC 
1.0A78E 00 
1.0173E 00 
9.9169E-01 
9.7157E-01 
9.5715E-01 
9.A823E-01 
9.AA30E-01 
9.AA63E-01 
9.A837E-01 
9.5A62E-01 
9.6250E-01 
9.7120E-01 
9.8001E-01 
9.8836E-01 
9.9581E-01 
1.0021E 00 
1.0070E 00 
1.0105E 00 
1.0127E 00 
1.0136E 00 
1.0135E 00 
1.0126E 00 
l.OlllE 00 
1.0092E 00 
1.007QE 00 
1.00A9E OO 
1.0028E 00 
1.00lOE OO 
9.99A8E-01 
9.9828E-01 
9.97A2E-01 
9.9689E-01 
9.9666E-01 
9.9669E-01 


OUT 

0.0 

5.8666E-03 
3.7258E-02 
1.0058E-01 
1.9206E-01 
3.0387E-01 
A.2726E-31 
5.5396E-01 
6.769AE-01 
7.9071E-01 
8.9133E-01 
9.763AE-01 
1.0AA6E OC 
1.0961E 00 
1.1317E 00 
1.1529E 30 
1.1617E 00 
1.1601E 00 
1.1505E 00 
1.1350E OO 
1.1158E 00 
1.0946E 00 
1.0729E 00 
1.0520E 00 
1.0328E 00 
1.0159E 00 
1.0018E 00 
9.9052E-01 
9.8215E-01 
9.76A8E-01 
9.732AE-01 
9.72C8E-01 
9.7262E-01 
9.7AA9E-01 
9.7729E-01 
9.8063E-01 
9.8A37E-01 
9.8810E-01 
9.9166E-01 
9.9A90E-01 
9.9773E-01 
I.OOOIE 00 
1.0019E 00 
1.0033E 00 
1.00A2E 00 
1.00ATE OO 
1.00A8E 00 
1.00A7E 00 
1.0043E 00 
1.0038E 00 
1.0032E 00 
1.0026E 00 
1.0019E 00 
1.0013E 00 
1.0008E 00 


OUT 

0.0 

A.6935E-03 

2.9823E-02 

8.0610E-02 

1.5A27E-01 

2. AA9AE-01 

3. A6C7E-01 
A.515AE-01 
5.5613E-01 
6.5572E-01 
7.A727E-01 
8.2873E-01 
9.9895E-01 
9.57A9E-01 
1.00A5E 00 
1.0A07E 00 
1.0668E OC 
1.08A2E 00 
1.09A1E 00 
1.0977E 00 
1.0963E 00 
1.0911E 00 
1.0832E 00 
1.0736E 00 
1.0630E 00 
1.0522E 00 
1.0A16E 00 
1.0316E 00 
1.0226E OC 
1.01A6E 00 
1.0C78E 30 
1.0023E 00 
9.9790E-01 
9.9A59E-01 
9.9227E-01 
9.9081E-01 
9.9007E-01 
9.8992E-01 
9.9C23E-01 
9.9089E-01 
9.9179E-01 
9.9283E-01 
9.939AE-01 
9.9506E-01 
9.9613E-01 
9.9712E-01 
9.9801E-01 
9.9878E-01 
9.99A3E-01 
9.9996E-01 
l.OOOAE 30 
1.0007E on 
1.0009E 30 
l.OOlOE 00 
l.OOlOE 00 


Computer output data for K = 55. 50, 45, and 36. 
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s.soSJe^oo 

5.6000E 00 
5.7000E 00 
5.8000E 00 
5.9000E 00 
6.0C00E 00 
6.1000E 00 
6.2000E 00 
6.3000E 00 
6.4000E 00 
6.5000E 00 
6.6000E 00 
6.700CE 00 
6.8000E 00 
6.9000E 00 
7.0000E 00 
7.1000E 00 
7.2000E 00 
7.300CE 00 
7.4000E 00 
7.5000E 00 
7.6OC0E 00 
7.7000E 00 
7.8000E 00 
7.9000E 00 
8.0000E 00 
8.1000E 00 
8.2000E 00 
8.3000E 00 
8.4000E 00 
8.5000E 00 
8.6000E 00 
8.7000E 00 
8.8C00E 00 
8.9C00E 00 
9.0000E 00 
9.1000E 00 
9.2000E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000E 00 
9.7000E 00 
9.8000E 00 
9,9000E 00 
l.OOOOE 01 


K=50 

9.983^E-01 
9.9823E-01 
9.9821E-01 
9.9829E-01 
9.9846E-01 
9.9868E-01 
9.9893E-01 
9,9920E-01 
9.9945E-01 
9.9969E-01 
9.9989E-01 
I.OOOIE 00 
1.0002E 00 
1.0003E 00 
1.0003E 00 
1.0004E 00 
1.0004E 00 
1.0003E 00 
1.0003E 00 
1.0003E 00 
1.0n02E 00 
1.0002E 00 
I.OOOIE 00 
I.OOOIE 00 
l.OOOOE 00 
9.9998E-01 
9.9995E-01 
9.9994E-31 
9.9992E-01 
9.9992E-01 
9.9992E-01 
9.9993E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9997E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 


K=55 

9.96§i’E-01 
9.9729E-01 
9.9776E-01 
9.9828E-01 
9.9881E-01 
9.9931E-01 
9.9975E-01 
I.OOOIE 00 
1.0004E 00 
1.0006E 00 
1.0008E 00 
1.0008E 00 
1.0008E 00 
1.0008E 00 
1.Q007E 00 
1.0005E 00 
1.0004E 00 
1.0003E 00 
1.0002E 00 
I.OOOIE 00 
9.9997E-01 
9.9989E-01 
9.9984E-01 
9.9981E-01 
9.9980E-01 
9.9980E-01 
9.9981E-01 
9.9984E-01 
9.9986E-01 
9.9989E-01 
9.9993E-01 
9.9996E-01 
9.9998E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 


K=45 

OUT 

1.0003E 00 

9.9992E-01 

9.9962E-01 

9.9940E-01 

9.9926E-01 

9.9919E-01 I 

9.9917E-01 

9.9920E-01 

9.9927E-01 

9.9936E-01 

9.9946E-01 

9.9957E-01 

9.9968E-0i 

9.9978E-01 

9.9988E-01 

9.9996E-01 

l.OOOOE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

I.OOOIE 00 

l.OOOOE 00 

l.OOOOE 00 

l.OOOOE 00 

l.OOOOE 00 

9.9999E-01 

9.9998E-01 

9.9998E-01 

9.9997E-01 

9.9997E-01 

9.9997E-01 

9.9997E-01 

9.9998E-01 

9.9998E-01 

9.9998E-01 

9.9999E-01 

9.9999E-01 

9.9999E-01 

9.9999E-01 

l.OOOOE 00 


Table 2.2. (Continued) 


OUT 

l.OOlOE 00 
l.OOlOE 00 
1.0009E 00 
1.0008E 00 
1.0007E 00 
1.0006E 00 
1.0005E 00 
1.0004E 00 
1.0003E 00 
1.0002E 00 
I.OOOIE 00 
l.OOOOE 00 
9.9998E-01 
9.9995E-01 
9.9992E-01 
9.9990E-01 
9.9989E-01 
9.9989E-01 
9.9989E-01 
9.9990E-01 
9.9991E-01 
9.9992E-01 
9.9993E-01 
9.9994E-01 
9.9995E-01 
9.9996E-01 
9.9997E-01 
9.9998E-01 
9.9999E-01 
9.9999E-01 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE on 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
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If the design criteria is to restrict the peak overshoot to 
+ 0.05 then (2.9) gives 


0.05 
0.366 


AK/K 


14% 


Figure 2.5 shows the results of the extreme values of K = 43 
and K. = 57 which give ’the tolerable - limits of overshoot. 


III. MULTILOOP CONTROL SYSTEM SENSITIVITY 


A. INTRODUCTION 

The Kokotovic [ 6 ] method of sensitivity analysis is 
expanded to investigate gain and time-contant variations in 
multiloop systems using a procedure similar to the one 
developed for the basic feedback control system. This proce¬ 
dure is developed using,-the block diagram of a regulated 
power generator which incl-udes an amplifier-exciter, a series 
RC compensating circuit in the feedback loop, and a load as 
shown in Fig. 3.1. 

B. GAIN-CONSTANT SENSITIVITY COEFFICIENTS 

Terms are defined as follows: 

KiGi(s) = Direct transfer function for the amplifier- 


exciter , 


K 2 G 2 (s) = Transfer function for the generator 


K 3 G 3 (s) = Transfer function for the load 


FH (s) 


Feedback .. 


or transfer function 


The basic equations are as follows: 



K G. {s j KzG , (s) K 3 G'. (s) 


I+K 1 G 1 is) FH(s)-^KiGi (s)K 2 G 2 (s)K 3 G 3 (s; 


(3.1a) 
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Ez 

r“ 


(s) 


1 

1+KiGi(s)FH(s)+K 1 G 1 (s)K 2 G 2 (s)K3G3(s) 


(3.1b) 


El 


1+KiGi(s)FH(s) 

1+KiGi (s)FH(s)+KiGi(s)K2G2(s)K3G3(s) 


{3.1c) 


Bi 

— (s) 


KiGi (s)FH(s) 

1+KiGi (s) FH (s)+KiGi (s)K 2 G 2 (s)K3G3 (s) 


(3.Id) 


As before, capital letters denote the Laplace transform 
of the function described, except for K and F which are gain 
constants, and T which is a time constant. 

To obtain the system sensitivity function apply (2.1) 
to (3.1a) to give 


S 


K, 


RK 1 G 1 K 2 G 2 K 3 G 3 _ 

(I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 )^ 

(3.2) 


_C_ 

I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 

This sensitivity function corresponds to the error signal, 

E 2 , in the original multiloop system but with the controlled 
output, C, as the input to the system. This can be seen by 
comparing (3.2) to (3.1b) with C replacing R. As in the 
basic control system of Section II, the sensitivity model for 
gain parameters can be realized by cascading two identical 
systems and using the error signals in the sensitivity model 
as the desired gain sensitivity coefficients as shown in 
Fig. 3.2. 


The load-gain sensitivity function is derived by applying 
(2.1) to (3.1a), yielding 
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Fl£.3.1. Multiloop system. 



























RK 1 G 1 K 2 G 2 K 3 G 3 (1+KiGiFH) 


^ (I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 ) 

C (1+KiGiFH) 
I+K 1 G 1 FH+K 1 G 1 K 2 G 2 K 3 G 3 


(3.3) 


This sensitivity function corresponds to the error signal, 
El, as can be seen by comparing (3.3) to (3.1c) with C 
replacing R as the input. 

For the feedbacJc-gain sensitivity function apply (2.1) 
to (3.1a), yielding 


RK 1 G 1 K 2 G 2 K 3 G 3 (-K 1 G 1 FH) 
I+K1G1FH+K1G1K2G2K3G3)^ 

C(-K 1 G 1 FH) 

1+K 1 G 1 FH+K 1 G 1 K 2 G 2 R 3 G 3 


(3.4) 


This sensitivity function corresponds to the negative of the 
feedback signal, Bi, as can be seen by comparing (3.4) to 
(3.Id) with C replacing R as before. 

By applying (2.1) to the system transfer function, 
sensitivity functions have been derived for the various gain 
parameters and the pickoff points for the sensitivity 
functions can be selected by inspection as shown in Fig. 3.2. 
No new components have been required and all sensitivity 
functions can be generated simultaneously. 


28 








29 


Fig.3.2. Multiloop sensitivity model. 






































c. 


TIME-CONSTANT SENSITIVITY COEFFICIENTS 


Since time-constant variations can also affect the 
response, sensitivity coefficients for LR and CR network, 
combinations can add helpful information to the analysis 
of a modeled system. For this derivation, let 

C _ KG _ KG(Ts+l) 

R 1+KGFH 1+Ts (KGF+1) ' 


where H = 
circuit. 


7 =—rr and T is the time constant of an RC feedback 
Ts+1 

Applying (2.1) to (3.5) yields 
-R(KG)^ FTs 

S = 

[1+Ts (KGF+1) ] ^ 


(3.6) 


C(-KGFTs) 1 

[1+Ts(KGF+1)] Ts+1 


and the feedback-gain sensitivity is 

^ _ C(-KGFTs) 

^F 1+Ts(KGF+1) 

Therefore, by combining (3.5) and (3.6) 

^T ^F ^fs+r^ 


(3.7) 


(3.8) 


This indicates that the RC time-constant sensitivity can 
be generated by using the feedback-gain sensitivity as the 
input to a block with elementary transfer function as shown 
in Fig. 3.2. 

In the direct path there are time-constants associated 
with the generator and load inductance and resistance 
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combinations. The effect of variations of these time- 
constants will add information to the analysis of the system. 
For this derivation, let 


C = kg ^ K 

R 1+KGFH 1+Ts+KFH 

where G = j -, and T is the time-constant of a direct path 
LR component. Applying (2.1) to (3.9) yields 


„ -RKTs 

(1+Ts+KFH)^ 

C(-Ts) 

(1+Ts+KFH) 

and the direct-gain sensitivity is 

^ C(Ts+l) 

K (1+Ts+KFH) 

Therefore, by combining (3.10) and (3.11) 


(3.10) 


(3.11) 


= ®K (fiS> 

As with the feedback RC time-constant sensitivity, only one 
elementary block has to be added in series with the direct- 
gain sensitivity model to complete the sensitivity model as 
shown in Fig. 3.2. 

D. APPLICATION 


One of the basic problem.s of preliminary design analysis 
is to estimate the influence of each design parameter on the 
system response and decide how improvement of the system 
can be assured. Sensitivity coefficients can be helpful xn 
solving this problem. To illustrate this procedure, a basic 
set of design parameters for a regulated generator is taken to 
be amplifier gain, exciter time-constant, load-gain. 
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1.5 



( 3EG ) 

F16.3.3. Multiloop system response and 
functions. 
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TIME ( SEG ) 

Fig,3.A. AiDplifler-exclter sensitivity functions. 
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TIIl£(bEC ) 


Fi6.3*5* liultiloop response for = 40 and = 60. 
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Pi6.3»6. Load and feedback gain sensitivity functions. 
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TIM£(Si!.C) 

Fig.3*7* Multiloop response for 1.1 and F= 2.2. 
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feedback-gain, and a compensator time-constant. Figure 3.3 
is a computer graph of the time response to a unit-step 
input and the five design-parameter sensitivity coefficients. 

At a glance, it is apparent from Fig. 3.3 that Ki and 

Ti have relatively little influence on the response. Table 

3.1 shows the maximum value of S„ is only 0.05 which is less 

^ 1 

than 10% the value of S„ . Figure 3.4 shows S„ and S^-, 

3 1 J-1 

where the scale has been increased by 10 to show the general 
shape of the curves. Figure 3.5 shows the time response for 
Kj = 40 and Kj =60. Even for this 50% increase, there is 
little change in the performance. Since the response is 
very "insensitive" to Kj and Tj, little improvement can be 
achieved by adjusting these parameters and their tolerances 
can be large in relation to the other parameters. It should 
be noted that the system is practically insensitive to small 
variations of resistance in the exciter because the time- 
constant and gain both either increase or decrease by the 
same amount. The sensitivities have the same relative 
amplitude but differ in sign and therefore the effect on the 
time response is canceled. 

In Fig. 3.6 the sensitivity coefficients for the load- 
gain and feedback-gain constants are plotted. Table 3.2 
shows that their magnitudes are almost equal but are 180 
degrees out of phase. Figure 3.7 is a graph of the time 
response with their original values and the response for a 
10% increase in both parameters. Although the response is 
sensitive to such changes, the plots are almost superimposed 
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due to the one being the negative of the other. This could 
be favorable equalization if aging or temperature cause 
these gain parameters to vary by the same percentage. 

Figure 3.3 shows that the feedback-gain and the compen¬ 
sator RC time-constant have almost the same sensitivity but 
the time-constant sensitivity function is delayed by one 
second. 

Since fractional parameter variations are considered, 
the permissible percentage of a parameter change can be 
determined from the allowable system response variation. 

It has been shown that any general form of a control 
system that can be transformed into a system of primary 
loops can be analyzed in the same fashion as developed in 
this section. The sensitivity to variations in direct- 
and feedback-gain constants are available as error or feed¬ 
back signals in the sensitivity model and time-constant 
sensitivities require addition of elementary transfer 
functions as shown in Fig. 3.2. 
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Tll'oi!. 


0.0 

l.OOOOE-01 
2.0000E-01 
3,0000E-01 
4.0000E-01 
5.0000E-01 
6,0000E-01 
T.OOOOE-Oi 
8.0000E-01 
9,0000E-0l 
l.OOOOE 00 
I.IOOOE 00 
1.2000E 00 
1.3000E 00 
1.4000E 00 
1.5000E 00 
1.6000E 00 
1.7000E 00 
1.8000E 00 
l,9000E 00 
2.0000E 00 
2.1000E 00 
2.2000E 00 
2.3000E 00 
2.4000E 00 
2*5000E 00 
2.6000E 00 
2.7000E 00 
2.8000E 00 
2.9000E 00 
3.0000E 00 
3.1000E 00 
3.2000E 00 
3.3000E 00 
3.4000E 00 
3.5000E 00 
3.6000E 00 
3.7000E 00 
3.8000E 00 
3.9000E 00 
4,0000E 00 
4.1000E 00 
4.2000E 00 
4.3000E 00 
4.4000E 00 
4,5000E 00 
4.6000E 00 
4.7000E 00 
4.8000E 00 
4*9000E 00 
5.0000E 00 
5.1000E 00 
5.2000E 00 
5.3000E 00 
5.4000E 00 


3Ki 


Sti 


0.0 

5.2083E-04 

2*9l89E-03 

7.4409E-03 

1.4014E-02 

2.2527E-02 

3.2866E-02 

4.4918E~02 

5.8570E-02 

7,37l3E-02 

9.0238E-02 

1.0804E-01 

1.2701E-01 

1.4705E-01 

l,6807E-01 

1.8995E-01 

2.1262E-01 

2.3597E-01 

2.5993E-01 

2.8440E-01 

3.0930E-01 

3.3455E-0X 

3.6007E-01 

3.8580E-01 

4,ll66E-0l 

4,3759E-0l 

4. 6351E-01 

4.8937E-01 

5.1512E-01 

5.4069E-01 

5.6603E-01 

5,9lllE-0l 

6*l586E-0l 

6.4025E-01 

6.6424E-01 

6,8779E-01 

7.1086E-01 

7,3344E-0l 

7.5548E-01 

7,7696E-0l 

7,9786E-0l 

8.1816E-0 

8.3784E-0 

8.5688E-0 

8.7528E-0 

8,930lE-0. 

9.1006E-01 

9,2644E-01 

9,4214E-0l 

9.5715E-01 

9,7147E-01 

9.8510E-01 

9,9804E-Ol 

1.0103E 00 

1.0219E 00 


? 


0*0 

3.2895E-04 
1.2793E-03 
2.4687E-03 
3.7582E-03 
5.1151E-03 
6,5283E-03 
7,9903E-03 
,4944E-03 
•1034E-02 
1.2603E-02 
1.4194E-02 
1.5801E-02 
1.7420E-02 
l,9043E-02 
2.0664E-02 
2.2280E-02 
2•3e83i-02 
2.5470E-02 
2.7035E-02 
2.8574E-02 
3.0082E-02 
3.1556E-02 
3.2991E-02 
3.4385E-02 
3.5733E-02 
3.7032E-02 
3.82 81E-02 
3.9476E-02 
4.0615E-02 
4,l696E-02 
4.2717E-02 
4.3677E-02 
4.4574E-02 
4.5408E-02 
4,6l78E-02 
4,6881E-02 
4,7520E-02 
4.8093E-02 
4,8600E-02 
4.9042E-02 
4.9419E-02 
4.9731E-02 
4.9980E-02 
5*0168E-02 
5.0293E-02 
5.0357E-02 
5,0362E-02 
5.0309E-02 
5.0203E-02 
5,0039E-02 
4.9825E-02 
4.9561E-02 
4.9247E~02 
4.8887E-02 


0,0 

-3.2650E-04 
-1.2578E-03 
-2,4017E-03 
-3.6l6iE-03 
'4.8670E-03 
-6.1426E-03 
-7.4349E-03 
-8.7368|-03 
-1.0042E-02 
-l,l343E-02 
-1.2635E-02 
-1.3910E-02 
-1.5165E-02 
-1.6393E-02 
-l,7590E-02 
-1.8751E-02 
-1.9872E-02 
-2.0948E-02 
-2.1977E-02 
-2.2954E-02 
-2.3877E-02 
-2.4742E-02 
•2.5549E-02 
-2*6294E-02 
-2.6976E-02 
-2,7594E-02 
-2.8146E-02 
•2.8631E-02 
-2*9048E-02 
*2,9399E-02 
-2.9681E-02 
•2.9896E-02 
•3,0044E-02 
-3.0126E-02 
■3.0142E-02 
-3,0092E-02 
-2,9980E-02 
-2.9805E-02 
-2.9570E-02 
•2.9276E-02 
'2,8925E-02 
-2.8520E-02 
•2,8062E-02 
-2.7553E-02 
-2.6997E-02 
-2.6393E-02 
-2. 5746E-02 
•2.5059E-02 
'2,4335E-02 
-2.3572E-02 
'2.2779E-02 
-2.1955E-02 
■2.1103E-02 
'2.0226E-02 


Table 3»la* Computer output data for C, 


3? 



TIM 

5.5066E00 
5.6000E 00 
5*7000E 00 
5*8000| 00 
5.9000E 00 
6.0000E 00 
6.1000E 00 
6.2000E 00 
6.3000E 00 
6.4000E 00 
6.5000E 00 
6.6000E 00 
6.7000E 00 
6*8000E 00 
•9000E 00 
00 
00 
00 
00 


I 


• OOOOE 

• lOOOf 


7.2000 
7.3000 
7.4000E 00 
7.5000E 00 
7.6000E 00 
7.7000E 00 
7.8000E 00 
7.9000E 00 
8.OOOOE 00 
8 .I 0 OOE 00 
8.2000E 00 
8.3000E 00 
8.4000E 00 
8.5000E 00 
8.6000E 00 
8.7000E 00 
8.8000E 00 
8.9000E 00 
9.OOOOE 00 
9.1000E 00 
9.2000E 00 
9.3000E 00 
9.4000E 00 
9.5000E 00 
9.6000E 00 
9.7000E 00 
9.8000E 00 
9.9000E 00 
1.OOOOE 0 
I.OIOOE 0 
1.0200E 0 
1.0300E 0 
1.040QE 0 
1.0500E 0 
i.0600E 0 
.07001 0 
1.0800E 0 
1.0900E 0 


.0328E 00 
.0430E 00 
.0526E 00 
•0615E 00 
.0698E 00 
.0775E 00 
.0845E 00 
.0910E 00 
.0969E 00 
.1021E 00 
.1069E 00 
.llllE 00 
.1148E 00 
.1179E 00 
.1206E 00 
.1228E 00 
.1246E 00 
.1259E 00 
.1268E 00 
.1273E 00 
.1274E 00 
.1271E 00 
•1266E 00 
.1256E 00 
.1244E 00 
.1229E 00 
.t21lE 00 
.1191E 00 
.I168E 00 
.I143E 00 
.1116E 00 
.1087E 00 
.1056E 00 
.1024E 00 
.0991E 00 
.0956E 00 
.0920E 00 
.0883E 00 
.0845E 00 
•0806E 00 
.0767E 00 
.0728E 00 
.0688E 00 
.0647E 00 
.0607E 00 
.0567E 00 
.0527E 00 
.0486E 00 
.0446E 00 
.0407E 00 
.0368E 00 
.0329E 00 
.0291E 00 
.0253E 00 
.0216E 00 


4.8484E-02 
4.8035E-02 
4.7550E-02 
4.7024E-02 
4.6467E-02 
4.5875E-02 
4.5251E-02 
4.4601E-02 
4.3923E-02 
4.3223E-02 
4.2501E-02 
4.1759E-02 
4.0999E-02 
4.0227E-02 
3.9443E-02 
3.8449E-02 
3. 7846E-02 
3.7040E-02 
3.6225E-02 
3.5413E-02 
3.4599E-02 
3.3787E-02 
3.2978E-02 
3.2178E-02 
3.1382E-02 
3.0596E-0| 
2.9821E-02 
2.9058E-02 
2.8305E-02 
2.7569E-02 
2.6848E-02 
2.6145E-02 
2.5457E-02 
2.4791E-02 
2.4144E-02 
2.3515E-02 
2.2910E-02 
2.2327E-02 
2.i766E-02 
2.1228E-02 
2.0713E-02 
2.0222E>02 
1.9754E-02 
1.9313E-02 
1.8897E-02 
1.8503E-02 
1.8134E-02 
1.7792E-02 
1.7472E-02 
1.7179E-02 
1.6908E>02 
I.6664E-02 
I.644jE-02 
L.6244E-0r 


®Ti 

-1.9329E-02 

-1.8408E-02 

-1.7475E-02 

-1.4597E-02 
-1.3620E-02 
-1.2641E-02 
-1.1660E-02 
-1.0681E-02 
-19.7037E-03 
-8.7310E-03 
-7.7655E-03 
-6.8106E-03 
-5.8688E-03 
-4.93921-03 
-4.0244E-03 
-5*12856-03 
-2.2473E-03 
-1.3897E-03 
-5.5103E-04 
2.6421E-04 
1.0564E-03 
1.8213E-03 
2.5623E-03 
3.2745E-03 
3.9590E-03 
4.6156E-03 
5.2448E-03 
5.8424E-03 
6.4102E-03 
6.9460E-03 
7.4535E-03 
7.9278E-03 
8.3706E-03 
8.7855E-03 
9.165fE-03 
9.51481-03 
9.8346E-03 

1.0610E-02 
1.0810E-02 
I.O979E-02 
1.1118E-02 
1.1233E-02 
1.1320E-02 
1.1378E-02 
1.1413E-02 
421E-02 
.407E-02 
366E-02 
,3051-0 
,22lE-0 


{: 


1.6069E-02 1.1116E-02 


Table 3.1a. (Continued) 
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~ui e. 
-01 7 

: 8 } J 


00 

00 

00 

00 

00 

00 


TIM£. 

0.0 

l.OOOOE-01 
2.0000E-01 
3.0000E-0 
4.0000E 
5.0000E 
6.0000E-01 
?• OOOOE-Ol 
8.0000E-01 
9.0000E-01 
l.OOOOE 00 
I.IOOOE 
1.2000E 
1.3000E 
1.4000E 
1.5000E 
1.6000E 
1.7000E 00 
1.8000E 00 
1.9000E 00 
2.0000E 00 
2.1000E 00 
2.2000E 00 
2.3000E 00 
2.4000E 00 
2.5000E 00 
2.6000E 00 
2.7000E 
2.8000E 
2.9000E 
3.0000E 
3.1000E 
3.2000E 
3.3000E 
3.4000E 
3.5000E 
3.6000E 
3.7000E 
3.8000E 
3.9000E 
4.0000E 
4.1000E 00 
4.2000E 00 
4.3000E 00 
4.4000E 00 
4,5000E 00 
4.6000E 00 
4.7000E 00 
4.8000E 00 
4.9000E 00 
5.0000E 00 
5.1000E 00 
5.2000E 00 
5.3000E 00 
5.4000E 00 


00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 


0*0 

5.2081E-04 
2.9185E-03 
7,4408E-03 
•4014E-02 
.2527E-02 
3.2866E-02 
4.4918E-02 
5.8570E-02 
7.3713E-02 
9,0238E-02 
1.0804E-01 
1.2701E-01 
l,4705E-01 
l.6807E-0i 
1.8995E-01 
2.1262E-01 
2.3597E-01 
2,5993E-0l 
2.8440E-01 
3,0930E-0l 
3.3455E-0i 
3.6007E-0i 
3.8580E-01 
4.1166E-01 
4,3759E-0l 
4.6351E-01 
4.8937E-01 
5.1512E-01 
5.4069E-01 
5,6603E-01 
5.9111E-01 
6.1586E-01 
6,4025E-0l 
6.6424E-01 
6.8779E-01 
7,1086E-0l 
7.3344E-01 
7.5548E-01 
7,7696E-0l 
7,9786E-0l 
8.1816E-01 
8.3784E-0I 
8.5688E-01 
8.7528E-01 
8.9301E-01 
9.1006E-01 
9,2644E-01 
9,42l4E-0l 
9.5715E-01 
9.7147E-01 
9,85i0E-0l 
9.9804E-01 
1,0103E 00 
1.0219E 00 


SKj 

- 0,0 

5,2079E-04 
2,9l77E-03 
7,4349E-03 

3.2718E-02 
4.4627E-02 
5,8053E-02 
7,2863E-02 
8,8919E-02 
l,0609E-0l 
1.2424E-0 
l,4323E-0 
l,6293E-0^ 
1.8323E-0I 
2.0398E-0 
2.2507E-0 
2,4638E-0_ 
2.6780E-01 
2.8921E-01 
3.1050E-01 
3,3158E-01 
3.5235E-01 
3.7271E-01 
3.9257E-0i 
4,1186E-01 
4,3051E-01 
4.4842E-01 
4.6555E-01 
4.8183E-01 
4,9721E-0l 
5.1162E-01 
5.2505E-01 
5.3743E-01 
5,4874E-01 
5.5895E-01 
5.6803E-01 
5,7597E-01 
5.8275E-01 
5,8837E-01 
5.928iE-0l 
5,9608E-0l 
5,9819E-0l 
5.9913E-01 
5.9892E-01 
5.9758E-01 
5.9513E-01 
5.9158E-01 
5.8696E-01 
5,8130E-0l 
5,7463E-01 
5.6698E-01 
5.5839E-01 
5,4890E-0l 


Sf 

u, u 

-1.9186E-04 
-1.6383E-03 
-4, 9661E-03 

:l:9illi:8l 

-2,6l90E-02 
-3.6636E-02 
-4.8559E-02 
-6.1829E-02 
-7,6317E-02 
-9,1894E-02 
-1.0843E-01 
-1.2581E-01 
-l,4389E-01 
-1.6256E-01 
-1.8170E-01 
-2.0119E-01 
-2.2091E-01 
-2.4076E-01 
-2.6063E-01 
-2.8042E-01 
-3.0002E-01 
-3,1935E-0l 
-3.3832E-01 
-3.5684E-0i 
-3,7483E-01 
'3,9222E-01 
•4.0895E-01 
-4.2494E-01 
-4,4014E-01 
-4.5449E-01 
-4.6795E-01 
-4,8047E-01 
-4.9202E-01 
•5.0256E-01 
-5.1206E-01 
•5.2051E-01 
•5.2788E-01 
•5.3415E-01 
■5.3933E-01 
-5.4339E-01 
-5.4635E-01 
•5,4821E-01 
-5,4896E-0l 
•5.4863E-01 
•5.4722E-01 
-5.4476E-01 
-5.4127E-01 
•5.3675E-01 
•5.3126E-01 
•5.2480E-0i 
•5.1742E-0i 
-5.0915E-01 
•5.0001E-01 


sth 

u,u 

-4,2725E-06 

-7,8408E-05 

-3.7232E-04 

-kSslii-SI 

-4,1039E-03 
-6,6987E-03 
-1.0113E-02 
-i,4404E-02 
-1.9608E-02 
-2.5750E-02 
-3.2838E-02 
-4,0866E-02 
-4,9819E-02 
-5,9671E-02 
-7.0385E-02 
-8.1918E-02 
-9.4222E-02 
-1.0724E-01 
-1.2091E-01 
-1.3516E-01 
-l,4993E-0l 
-1,6516E-01 
-1.8075E-01 
-1.9664E-01 
-2.1276E-01 
-2.2903E-01 
-2,4538E-01 
-2,6l72E-0l 
-2,7800E-01 
-2.9413E-01 
-3.1005E-01 
-3.2569E-01 
-3.4098E-01 
-3.5587E-01 
-3.7030E-01 
-3.8421E-01 
-3.9754E-01 
-4.1026E-01 
-4,2231E-01 
-4.3365E-01 
-4.4425E-01 
-4.5406E-01 
-4.6307E-01 
-4.7123E-01 
-4.7854E-01 
-4,8496E-0l 
-4.9049E-01 
-4.9511E-01 
-4,9882E-01 
-5.0160E-01 
-5.0346E-01 
-5.0439E-01 
-5.0441E-01 


Table 3•lb. Computer output data for C, , Sp, and S^g. 
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00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

01 

01 

01 

01 

01 

01 

01 

01 

01 

01 

) . 


•0328e 00 
•0430E 00 
.0526E 
•0615E 


00 

00 


•0698E 00 
.0775E 00 
•0845E 00 
•0910E 00 
•0969E 00 
•1021E 00 
•1069E 00 
•llllE 00 
•1148E 00 
.1179E 00 
•1206E 00 
•1228E 00 
•1246E 00 
.1259E 00 
•1268E 00 
•1273E 00 
.1274E 00 
.1271E 00 
'166E 00 
»56E 00 
l44E 00 
•1229E 00 
00 
00 

•1168E 00 
•1143E 
•1116E 


:H4li 


00 

00 


•1087E 00 


•1056E 
•1024E 
.0991E 
.09561 
•0920E 
.0883E 
•0845E 00 
•0806E 00 
.0767E 00 
.0728E 00 
•0688E 00 
•0647E 00 
•0607E 00 
.0567E 00 
.0527E 
•0486E 
.0446E 
.0407E 00 
•0368E 00 
•0329E 00 
•0291E 00 
.0253E 00 
•0216E 00 


00 

00 

00 

00 

00 

00 


00 

00 

00 


-ui - 


SK3 

5.3854E-01 
5.2736E-01 
5.1539E-01 
5.0269E-01 
4.8929E-01 
4.7525E-01 
4.6060E-01 
4.4540E-01 
4.2969E-01 
4. 1352E-0i 
3.9694E-01 
3,7999E-01 
3.6273E-01 
3.4520E-01 
3.2744E-01 
3.0951E-01 
2.9144E-01 
2.7329E-01 
2.5509E-01 
2.3690E-01 
2.1874E-01 
2.0066E-01 
1.8270E-0 
i.6490E 
1.4728E 
1.2990E-01 
1,1277E-01 - 
9.5940E-02 
7.9426E-02 
6.3260E-02 
4,7470E-02 
3.2081E-02 
1.7114E-02 
2.5921E-03 
-1.1465E-02 
-2.5040E-02 
-3.8116E-02 
-5.0678E-02 
-6.2716E^02 
-7.4216E-02 
-8.5169E-02 
-9.5567E-02 
-1.0540E-01 
-1.1468E-01 
-1.2338E-01 
-1.3151E-01 
-1.3907E-01 
-1.4606E-01 
-1.5247E-01 
-1.5833E-01 
-1.6362E-01 
-1.6836E-01 
-1.7255E-01 
-1.7620E-01 
-i.7932E-01 


Sp 

4.9006E-01 

4,7932E-0 

4.6784E-0 

4.5566E-0: 

4,4282E-0 

4,2937E-0: 

4,1535E-0 

4.0079E-0 

3.8576E-0 

3,7029E-0 

3,5443E-0 

3.3823E-0 

3,2173E-0 

3,0497E-0 

2.8800E-0 

2.7086E-0 

2.536QE-0 

2.3625E-0 

2.1887E-0 

2.0148E-0 

1.8414E-0 

l,6687E-0 

1.49721-0, 

1.32721-0 . 

1.1590E-0 

9.9302E-02 

8.2952E-02 

6.6882E-02 

5.1121E-02 

3.5691E-02 

2.0622E-02 

5.9357E-03 

8.3437E-03 

2.2199E-02 

3.5609E-02 

4.8554E-02 

6.1026E-02 

7.3006E-02 

8.4481E-02 

9.5444E-02 

1.0588E-01 

1.1579E-01 

1.2516E-01 

1.3399E-01 

1.4227E-01 

1.5001E-01 

1.5720E-01 

1.6385E-01 

1.6995E-01 

1.7551E-01 

1.8053E-01 

1.8502E-01 

1.8899E-01 

1.9245E-01 

1.9539E-01 


Sth 

5.0352E-01 

5.0172E-01 

4.9904E-01 

4.9549E-01 

4.9108E-01 

4.8584E-01 

4.7980E-01 

4.7296E-01 

4.6537E-01 

4.5705E-0 

4.4803E-0 

4.3834E-0 

4.2802E-0 

4.1710E-0 

4.0561E-0 

3.9359E-01 

■3.8107E-01 

3.6810E-01 

3.5471E-01 

3.4094E-01 

3.2683E-01 

3.1242E-01 

2.9774E-01 

2.8283E-01 

2.6773E-01 

2.5247E-01 

2.371QE-01 

2.2165E-0I 

2.0616E-01 

1.9066E-01 

1.7518E-01 

1.5976E-01 

1.4442E-01 

1.2921E-01 

1.1415E-01 

9.9269E-02 

8.4595E-02 

7.0153E-02 

5.5971E-02 

4.2070E-02 

2.8474E-02 

1.5205E-02 

2.2815E-03 

1.0278E-02 

2.2456E-02 

3.4237E-02 

4.5607E-02 

5.6552E-02 

.7062E-02 

.7126E-02 

8.6736E-02 

9.5883E-02 

1.0456E-01 

1.1277E-01 

1.2050E-01 




intlnued ) 
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IV. EXAMPLES OF DIFFERENTIAL EQUATION SENSITIVITY ANALYSIS 


A. INTRODUCTION 

A general method of determining the sensitivity functions 
for the coefficients of a second-order differential equation 
is derived for a mass, spring, and damper system. The 
procedure is then applied to the study of a sixth-order 
differential equation, the solution of which is initially 
unstable. The coefficients are then adjusted to yield a 
stable solution using the results of sensitivity analysis. 

b'. second-order differential equation sensitivity to 

COEFFICIENT VARIATIONS 

The elementary mass-spring-damper system shown in Fig. 

4.1 is described by the following differential equation 

MX + BX + KX = u(t) (4.1) 

Taking the Laplace transform and solving for X(s) gives 

X(s) = - (4.2) 

Ms +BS+K 

For convenience, let M = 1. Applying (2.1) to (4.2) yields 

S = ^ = X(s) -) (4.3) 

(s +BS+K) s +BS+K 

By comparing (4.3) to (4.2) it can be seen that the sensi¬ 
tivity model is a replica of the original system with the 
controlled output as the input to the model. The sensitivity 
coefficient is measured at the point corresponding to the 
feedback signal. 
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Fig.4.1. Blocli dlagran? of TDass-spring-damper syateir. 



-S 


K 


-Sb 


Fig.4.2. Second-order differential equation sensitivity 
model. 
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similarly, 


S = _ B s U(s) ^ (—HBs j 

(s^+Bs+K)^ s^+Bs+K 

Figure 4.2 is a block diagram of the complete sensitivity 
model which shows clearly that Sj^ and Sg are comparable to 
the feedback gain sensitivities as derived in Section II 
and III. 

The solution of (4.1) for a unit-impulse input is 

T _ B __ 

x(t) = - e 2 sin B^ t (4.5) 

/ - v 

^ J- 


Figure 4.4 is the graph of x(t) and the sensitivity coeffi¬ 
cients Sj^ and Sg for B = 0.6 and K = 0.6. The plot shows 
that the system is more sensitive to changes in K than B 
as verified by (4.5) . Table 4.1 is the computer output data 
for X/ and Sg. 

C. SIXTH-ORDER DIFFERENTIAL EQUATION SENSITIVITY TO COEFFI¬ 
CIENT VARIATIONS 


Consider the differential equation 


d®x ,d^x d^x d^x d^x ^x 

^ Id o d — 0 - ^ ^ ^ g ^ 


(4.6) 

u(t) 


dt‘ 


dt = 


dt‘ 


dt^ 


dt‘ 


dt 


where nominally a = 1, b = 3.4, c = 11.4, d = 9.4, e = 7.4, 
f = 0.575, g = 0.452. The solution when u(t) is a unit-step 
can be determined by DSL using only one block. For the 
sensitivity model the equation must be integrated in six 
steps to isolate eaph sensitivity coefficient. The procedure 
used for the second-degree equation can be expanded to give 
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^(s) 
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Fig.A.3. Sixth-order differential equation sensitivity model. 










































TIM (SEC ) 

Fig.4.4. Second-order systaiu response and sensitivity 
functions. 
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Fig.4.5. Sixth-order unstable system response and 
sensitivity functions. 

(a) . Response and sensitivity functions for 

d and f. 

(b) . Sensitivity functions for b and o. 

(c) . Sensitivity functions for e and g. 



Fig.4.5(a). X, Sjj , and Sf. 
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Flg.4.5(b). S-b and Sq . 



'Fig.4.5(c)* Sq and Sg. 
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Fig.4.6. Sixth-order stable system response and 
sensitivity functions. 

(a) . Response and sensitivity functions for 

d and f. 

(b) . Sensitivity functions for b and c. 

(c) . Sensitivity functions for e and g. 



miE:(SRC) 


Fig.4.6 (a). X, S< 3 , and Sf. 
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- 1.0 

Fig.4. 
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-2 


■ TIIffiOEC) 

Fl6.4.6(c). Sq and Sg. 


■TII4E(SEC) 
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a model similar to Fig. 4.2. The expanded sensitivity model 
is shown in Fig. 4.3. The solution of (4,6) is unstable as 
can be seen in Fig. 4.5(a). The sensitivity functions 
graphed in Fig. 4.5(b) indicate that the inner-loop feed¬ 
back gain parameters b and c influence the output relatively 
little. Going from the coefficient of the highest order 
term to the lowest, the influence increases but it is ob¬ 
tained with an increasing time delay. 

The problem of stabilizing the solution may be solved 
using sensitivity analysis by noting that is out-of-phase 
with X. By increasing the coefficient f, the response 
becomes stable as shown in Fig. 4.6(a) for f = 1.8. Any 
further adjustment of the performance can be achieved by 
referring to the sensitivity functions for the stabilized 
system as shown in Fig. 4.6. 
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TIMc 

0.0 

1 . 0000^-01 
2 . 0000 F -01 
3.0000F-01 
4.0000F-01 
5.00006-01 
6.00006-01 
7.00006-01 
8.00006-01 
9.00006-01 
1.00006 00 

1.10006 00 

1.20006 00 

1.30006 00 

1.40006 00 

1.50006 00 

1.60006 00 

1.70006 00 
1.B0006 00 

1.90006 00 
2.00006 00 
2.1000P 00 

2.20006 00 

2.30006 00 

2.40006 00 

2.50006 00 

2.60006 00 

2.70006 00 

2.80006 00 

2.90006 00 
3.00006 00 

3.10006 00 

3.20006 00 

3.30006 00 

3.40006 00 

3.50006 00 

3.60006 00 

3.70006 00 

3.80006 00 

3.90006 00 
4.00006 00 

4.10006 00 

4.20006 00 

4.30006 00 

4.40006 00 

4.50006 00 
4.60C06 00 

4.70006 00 

4.80006 00 
4.90C06 00 
5.0000'^ 00 

5.10006 00 

5.20006 00 

5.30006 00 

5.40006 00 


X 

0.0 

4.89906-03 
1.91856-02 
4.22296-02 
7.33866-02 
1.12006-01 
1.57426-01 
2.08986-01 
2.66016-01 
3.27886-01 
3.93946-01 
4.63566-01 
5.36136-01 
6.11066-01 
6.87766-01 
7.65706-01 
8.44346-01 
9.23196-01 
1.00186 00 
1.07976 00 
1.15646 00 
1.2317F 00 
1.30516 00 
1.37636 00 
1.44516 00 
1.51126 00 
1.57436 00 
1.63436 00 
1.69096 00 
1.74426 00 
1.79396 00 
1.83996 00 
1.8823E 00 
1.92096 00 
1.95586 00 
1.98706 00 
2.01456 00 
2.03846 00 
2.05876 00 
2.0756E 00 
2.08916 00 
2.09946 00 
2.10656 CO 
2.11076 00 
2.11206 CO 
2.11076 00 
2.10686 00 
2.10066 00 
2.09236 00 
2.08206 00 
2.O690E 00 
2.05606 00 
2.04086 00 
2.02426 00 
2.00666 00 


S6NSB 

- 0.0 

-9.69956-05 
-7.51876-04 
-2.45636-03 
-5.63C1P-C3 
-1.0622E-02 
-1.77136-02 
-2.71146-02 
-3.89766-C2 
-5.33856-02 
-7.03726-02 
-8.99126-02 
-1.11936-01 
-1.3632E-C1 
-1.629G«=-01 
-1.91506-01 
-2.21076-01 
-2.5377E-01 
-2.86026-01 
-3.21026-01 
-3.55776-01 
-3.90856-01 
-4.25946-01 
-4.60726-01 
-4.94876-01 
-5.28076-01 
-5.60036-01 
-5.90456-01 
-6.19056-01 
-6.45580-01 
-6.69816-01 
-6.91516-01 
-7.10496-01 
-3.26596-01 
-7.39666-01 
-7.49596-01 
-7.56296-01 
-7.59686-01 
-7.59746-01 
-7.56466-01 
-7.49856-01 
-7.39956-01 
-7.26826-01 
-7.10556-01 
-6.91256-01 
-6.69066-01 
-6.44106-01 
-6.16566-01 
-5.86616-01 
-5.54456-01 
-5.20276-01 
-4.84306-01 
-4.46766-01 
-4.07886-01 
-3.67896-01 


SENSK 

- 0.0 

-2.44016-06 
-3.80796-05 
-1.87896-04 
-5.78406-04 
-1.37466-03 
-2.77296-03 
-4.99436-03 
-8.27786-03 
-1.28746-02 
-1.90416-02 
-2.70346-02 
-3.71066-02 
-4.95006-02 
-6.44436-02 
-8.21486-02 
-1.02806-01 
-1.26576-01 
-1.53606-01 
-1.83996-01 
-2.17836-01 
-2.55166-01 
-2.96006-01 
-3.40346-01 
-3.88126-01 
-4.39286-01 
-4.93706-01 
-5.51236-01 
-6.11736-01 
-6.74986-01 
-7.40776-01 
-8.08856-01 
-8.78986-01 
-9.50866-01 
-1.02426 00 
-1.09876 00 
-1.17406 00 
-1.24986 00 
-1.32586 00 
-1.40176 00 
-1.47706 00 
-1.55156 00 
-1.62496 00 
-1.69686 00 
-1.76696 00 
-1.83496 00 
-1.90066 00 
-1.96376 00 
-2.02396 00 
-2.08096 00 
-2.13476 00 
-2.18496 00 
-2.23156 00 
-2.27426 00 
-2.31306 00 


Table 4.1. Computer output data for X, Sb» and 
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TIMI 


X 


S5NSR S5NSK 

•5 .5000*^ 

OC 

1.98805 

CC 

-3.27C35-C1 -2.34785 

00 

5.60C0'^ 

00 

1.9686E 

CC 

-2.85555-Cl -2.3784C 

00 

5.7000E 

00 

1.94855 

CC 

-2.4366E-C1 -2.40495 

00 

5.80C0P 

00 

1.928C5 

CC 

-2.01625-C1 -2.42715 

00 

5.90C0P 

OC 

1.9C725 

CC 

-1.59655-01 -2.44525 

00 

6.0000P 

OC 

1.88625 

OC 

-1.17985-01 -2.45915 

00 

6.1000P 

00 

1.06515 

CO 

-7.68155-02 -2.46885 

00 

6.200CF 

OC 

1.84415 

00 

-3.63825-02 -2.47455 

00 

6.300CE 

OC 

1.82335 

OC 

3.12685-C3 -2.47615 

00 

6.40005 

OC 

1.80285 

CC 

4.152C5-C2 -2.47395 

00 

6.50005 

OC 

1.78275 

CC 

7.86165-02 -2.46795 

00 

6.600CF 

00 

1.76315 

CC 

1.14255-ci -2.45825 

00 

6.7000F 

OC 

1.74415 

CC 

1.48275-Cl -2.44515 

00 

6.80005 

00 

1.72575 

00 

1.8C54F-C1 -2.42865 

00 

6.9000F 

OC 

1.70805 

CC 

2.10925-Cl -2.40905 

00 

7.00C05 

OC 

1.69125 

OC 

2.3931F-C1 -2.38655 

00 

7.1OC0F 

OC 

1.67525 

CO 

2.65615-01 -2.36125 

00 

7.2000E 

OC 

1.660CE 

CC 

2.09745-01 -2.33345 

00 

7.3000E 

OC 

1.64585 

CC 

3.11635-01 -2.30335 

00 

7.40CCF 

OC 

1.63265 

CC 

3.3123E-C1 -2.27125 

00 

7.5000E 

00 

1.62035 

00 

3.4P515-C1 -2.23725 

00 

7.60CCE 

OC 

1.60895 

CC 

3.6344F-C1 -2.20165 

00 

7.7000S 

00 

1.59865 

CO 

3.76C1F-C1 -2.16465 

00 

7.80005 

00 

1.58935 

CC 

3.86235-01 -2.12645 

00 

7.90C0F 

OC 

1.58105 

CC 

3.9412F-C1 -2.08745 

00 

8.0000E 

00 

1.57365 

CO 

3.9972F-C1 -2.04775 

00 

8.10C05 

OC 

1.56735 

CC 

4.C3C65-C1 -2.00755 

00 

8.20C05 

00 

1.56185 

OC 

4.Ci20F-Cl -1.96715 

00 

8.30CCP 

OC 

1.55745 

CO 

4.03225-Cl -1.92685 

00 

8.4000E 

00 

1.55385 

OC 

4.0C17E-01 -1.88665 

00 

8.50005 

OC 

1.551CE 

CC 

3.95165-01 -1.84685 

00 

8.60COE 

00 

1.54915 

CC 

3.88285-01 -1.80765 

00 

8.70C05 

00 

1.54805 

CC 

3.79625-01 -1.76925 

00 

3.80CC5 

00 

1.54775 

CC 

3.6930E-C1 -1.7317E 

00 

8.900CC 

00 

1.54805 

OC 

3.5743P-01 -1.69545 

00 

9.00005 

OC 

1.54915 

CC 

3.4413F-C1 -1.66035 

00 

«.10005 

00 

1.55C7E 

CC 

3.2952E-C1 -1.62665 

00 

9.20CC5 

OC 

1.55295 

OC 

3.13735-Cl -1.59445 

00 

9.3000E 

00 

1.55575 

CC 

2.9689F-C1 -1.56395 

00 

9.40005 

00 

1.55905 

CO 

2.7913E-C1 -1.53515 

00 

9.50CC5 

00 

1.56275 

CC 

2.6C57F-C1 -1.50815 

00 

9.60005 

00 

1.56675 

OC 

2.4130F-C1 -1.48305 

00 

O.700CE 

OC 

1.57125 

CC 

2.216CF-C1 -1.45985 

00 

9.8000'^ 

00 

1.57595 

CC 

2.C144F-01 -1.43875 

00 

9.90CC5 

00 

1.58095 

CC 

1.810C9-C1 -1.41955 

00 

1 .00005 

01 

1.58605 

OC 

1.6C4C5-C1 -1.40259 

00 

I.OIOOE 

01 

1.5*5145 

CO 

1.39765-Cl -1.38755 

00 

1.02C05 

01 

1.59695 

CO 

1.192C5-C1 -1.37455 

00 

1.03005 

01 

1.60245 

CO 

9.8811F-02 -1.36365 

00 

1.0400E 

01 

1.6081F 

CC 

7.87135-02 -1.35475 

00 

1.05005 

01 

1.6137F 

OC 

5.8999F-02 -1.34795 

00 

1.06005 

01 

1.61935 

00 

3.97635-02 -1.34295 

00 

1.0700E 

01 

1.62495 

OC 

2.1C905-C2 -1.3399E 

00 

1.08005 

01 

1.63035 

OC 

3.06125-03 -1.33875 

00 

1.0900E 

01 

1.63575 

CC 

-1.42525-02 -1.33925 

00 


Table 4.1. (Continued) 
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V. NONLINEAR SYSTEM SENSITIVITY 


A. INTRODUCTION 

While it is true that no physical systems are exactly 
linear, often adequate linear approximations can be used. 

In these cases, the Kokotovic [6] method of sensitivity 
analysis as developed in the preceding sections of this 
thesis can be applied. There are many situations, however, 
for which linear representations are invalid and nonlinear 
characteristics must be considered. 

Figure 5.1 is a block diagram of a nonlinear system 
containing a linear element represented by a transfer function 

K , where D = d/dt is the differential operator, and a 

nonlinear element, N, whose output is denoted by f(e). D is 
used instead of s in the linear transfer function because the 
Laplace transfomv and its inverse are generally not suitable 
for nonlinear analysis. 

The sensitivity to gain-constant variations is developed 
for the general nonlinear system and then applied to a 
specific system containing a saturation element as the non¬ 
linear component. 

B. SENSITIVITY COEFFICIENT FOR A GAIN FACTOR 

In order to designate the dependence of the variables in 
Fig. 5.1 the following notation is used: 
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Fis.5»l* Block diagrani of nonlinear system. 



Fig. 5 .2. Sensitivity model of nonlinear system. 
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r 


r(t) 


e 

= 

e(t,K) = r(t) - c(t,K) 


f (e) 

= 

output of nonlinear element N 


c 

= 

c(t,K) = r(t) - e(t,K) 


X 

= 

x(t,K) 


P(D) 

= 

,n ,q 

(5.1a) 

Q(D) 

= 

‘’n dt Vl ^^n-1 + ••• + + b. 

■ (5.1b) 


From the definition of gain-constant sensitivity coeffi¬ 
cient (2.1) 


_ 3c _ 3e 

K aK/K 3K 


(5.2) 


Using the differential equation for the linear element of the 
system gives 


KP(D) f(e(t,K)) = Q(D) c(t,K) = Q(D) [r (t)-e (t ,K) ] (5.4) 


Taking the logarithmic derivative, 3/3 in K = K 3/3 K, of 


(5.4) and 
KP (D) 


substituting (5.2) yields 
f(e(t,K)} + K^P (D) ) - = 

-KQ(D)|fi^ 


(5.5) 


Rearranging and simplifying (5.5) gives 


S 


K 


i® = KP (D) . K^P(D) 3f (e) 

3K Q(D) ^ Q(D) 3K 


(5.6) 


Figure 5.2 is a block diagram of the proposed sensitivity 
model for the nonlinear system. It will be shown that the 
describing equation of this model is identical to (5.6) . 
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From Fig. 5.2 


“ 9k/K " 

Substituting (5.2) in (5.7) yields 


(5.7) 


(5.8) 


Applying the transfer function for the linear element of 
the model gives 


^ ^ ^ Q(D)c + KQ(D) ||- 


-K'^P(D) ^ 


)e 3K -V.-/ gj^ 

Substituting (5.4) in (5.9) for Q(D)c and solving for S 
as defined by (5.2) gives 


(5.9) 


K 


'K 


K3e _ ,, P(D) . .. 2 P(D) 9f(e) 

~Tk - ^ qW ^ qTdT 3K 


(5.10) 


Equations (5.10) and (5.6) are the same and, hence, the model 
generates the gain-constant sensitivity coefficient of inter¬ 
est. The sensitivity model block represented by ^g ~ ^ '~ is th® 
instantaneous slope of the nonlinear function f(e). In many 
nonlinear systems a piecewise-linear approximation can be 
utilized. The sensitivity model will be the same as that of 
a linear system except that the gain changes as a function of 
the error signal in the original nonlinear system. An excimple 
for a piece wise-linear saturable element will make this more 
meaningful. 


C. APPLICATION 

The nonlinear component used as an example will be a 
limiter with a saturation characteristic as shown in Fig. 5.3. 
The linear plant used in Section II is used for comparison 
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of the time responses and gain sensitivity coefficients. 

The transfer function in differential operator notation 
becomes 

KP (p) ^ _50_ 

+ lOD^ + 24D 

The gain-constant sensitivity model as shown in Fig. 5.2 
requires a piecewise-linear function described by Fig. 5.4. 
It can be seen that is the slope of the nonlinear 

function f(e) and may be handled by logic statements in the 
DSL computer program as given in Appendix D. 

Figure 5.5 shows the time response and the sensitivity 
function for a unit-step input with a limited amplitude of 
error signal set at 0.1. The nonlinear element is saturat¬ 
ed until 4.8 seconds. During this period of saturation the 
gain representing the nonlinear element in the sensitivity 
model is zero, and therefore, the sensitivity coefficient 
and the output of the nonlinear system are one and the same 
as shown in Fig. 5.5 for the first 4.8 seconds. Figure 5.6 
shows the response and sensitivity for the limited amplitude 
of pi = -1.0 and pa = 1.0. This lets the system perform in 
a completely linear fashion and the effect of the nonlinear 
element may be realized by comparison with Fig. 5.5. The 
nonlinear response is slowed down, the maximum sensitivity 
to gain is greater but delayed by 3.8 seconds; also, the 
peak overshoot is greatly reduced. 

To verify the sensitivity coefficients for gain vari¬ 
ations, check points were chosen in the saturation region, 
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f*(e)=Pi, 
f(e)=P2, 
f(e)=e , 


?’i6*5«3* Characteristics of nonlinear N. 



9f(e) = 1 
3 e 

3 f (e ) = 0 
d e 

(e) = 0 
a e 


Fi6*5*^* Characteristics of sensitivity model 


e< Pi 
e>P2 

Pl^e < P2 


, Pl^e < P2 

, e <pi 

, e> P 2 


nonlinearity. 
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1.5 



-0.5 -^^^- 

Fi6«5.5. Response and sensitivity for nonlinear system 



Fig.5*6. Response and sensitivity for unlimited linear 


system. 


i- 
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peak-positive, and peak-negative portions of the sensitivity 
curve as shown in Fig. 5.7 for pi = -0.2 and pz = 0.2. The 
approximate change in the response due to a fractional change 
of 0.1 in the gain parameter, 

AC AK/K 

was calculated and.compared to the computer output for K = 55 
and K = 45 as given in the Table 5.2. The maximum error 
found was 0.5% for point 2 as presented in Table 5.1. 

Figure 5.8 shows the response for 14% variations in 
gain as was required for the linear system to constrain the 
peak overshoot to within + 0.05. The peak overshoot for 
the nonlinear system was constrained to + 0.014 for the same 
gain variations. This can be predicted from Fig. 5.8 since 

C. 

at the time of peak overshoot, the sensitivity function is 
much nearer zero than the sensitivity function for the linear 
system shown in Fig. 5.7. 
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1.5 



TIl«iE{SjtG) 

FIs. 5.7. Reaponae and aenaitivlty for nonlinear ayatem 
with P 3 _=*- 0.2 and P 2 = 0.2. 


Pt. 

T 

^50 

Sk 

C+AG 

G 55 

C-AG 

G 45 

1 

1.5 

.45215 

.45215 

.49737 

.49736 

.40694 

.40693 

2 

2.5 

.86763 

.84752 

.95238 

.94755 

.78288 

.78127 

3 

4.0 

1.0177 

-.21111 

.99659 

.99673 

1.0388 

1.0351 


Table 5.1. Comparison of actual outputs and outputs 
predicted from aenaitivlty coefficients. 
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8 . Gain tolerance of 14^ for overshoot of i 0.014. 
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K=50 


K=55 


K=45 


T IMF 

0.0 

l.OOOOE-Ol 
2.0000P-01 
3.0000P-01 
<*.00C0E-01 
•^.OOOOE-Ol 
6.0000E-01 
7,OOOOE-Ol 
8.0000E-01 
9.0000E-01 
l.OOOOE 00 
I.IOOOF 00 
1.2000E OC 
1.3000E OC 

l.AOOO«= 00 
1.5000E 00 
1.6000E 00 
1.7000E 00 

1. fiOOOE 00 
1.9000F OC 
2.0000E 00 
2.1000E 00 
2.2000F 00 
2.3000E CC 

2. AO00E OC 
2.500CE CC 
2.6000E 
2.7000= 
2.8000E 
2.9000E OC 
3.0000E 00 
3.100CE 00 
3.2000E 00 
3.3000E OC 

3. AOOOE CO 
3.500CE 00 
3.600CE 00 
3.7000E 00 
3.8000F 00 
3.9000E 00 
A.OOCOE 00 
A.IOOOE OC 
4.2000E 00 
A.3000E 00 
4.4000E OC 
A.'iOOOE 00 
A.60COE 00 
A.7000E 00 
A.8000P 00 
A.9000F 00 
•5.00COF 00 
5.10COE OC 
5.2OC0E 00 
5.3000E 00 
5.A0C0E 00 


CO 

OC 

OC 


SENSK 

C.O 

1.30A1E-C3 

8.302AE-C3 

2.2552E-02 

A.35A7E-02 

7.0099E-02 

1.009AE-01 

1.3A98E-C1 

1.7132E-01 

2.0930E-C1 

2.A8A3E-C1 

2.8837E-G1 

3.2886E-C1 

3.6972E-C1 

A.1085E-C1 

A.5215E-01 

^.9357E-C1 

5.3507E-C1 

5.7662E-01 

6.1821E-CI 

6.5983E-C1 

7.01A6E-C1 

7. A311E-C1 
7.8A76E-C1 
8.2A9AE-01 

8. A752E-C1 
0.39OfiE-Ol 
7.9771=-C1 
7.2738E-C1 
6.3A57E-C1 
5.2669e-Cl 
A.1113E-01 
2.9A68E-01 
1.8319F-01 
‘>.13A^E-C2 

-7.3A81E-C3 

-8.0616E-C2 

-1.3732E-C1 

-1.7736E-C1 

-2.01A7E-C1 

-2.1111E-C1 

-2.0823E-C1 

-1.95C9E-C1 

-1.7A12E-01 

-1.A77AE-C1 

-1.1826E-C1 

-e.7735E-02 

-E.7953E-02 

-3.03A5E-C2 

-5.9862E-03 

1.AA02E-02 

3.0A35E-C2 

A.2025E-C2 

A.9333E-C2 

5.2723E-02 


niT 

0.0 

1.3CA1E-C3 

8.302AF-03 

2.2552F-C2 

A.35A7F-C2 

7.0C99E-02 

1.009AF-C1 

1.3A98E-C1 

1.7132E-01 

2.093C=-C1 

2.A8A3F-01 

2.8837E-C1 

3.2886F-01 

3.6972=-Cl 

A.1085F-01 

A.5215E-C1 

A.9357F-C1 

5.35C7F-C1 

5.7662F-01 

6.1821F-C1 

6.5983F-C1 

7.01A6F-C1 

7.A311E-C1 

7.8A76F-C1 

8.26A1F-C1 

8.6763E-C1 

9.072CC-C1 

9.A366F-C1 

F.7577F-C1 

1.CC27F CO 

1.C239F CO 

1.0395E CC 

1.0A96F CO 

1.C5APF CC 

1.05F8F CO 

1.C533F CO 

1.CA83E CC 

I.CA16F CO 

1.C33PE 00 

1.0256F CO 

1.C177F CO 

1.C1C3F 00 

1.CC3PF 00 

9.98A5F-C1; 

9.9A25E-C1' 

9.912AF-C1 

9.893AF-C1 

9.88AAF-01i 

9.8839F-C1 

9.89C1F-01 

9.901AE-01 

9.F161F-C1 

9.9327E-C1 

9.9A97F-01 

9.9662F-01 


017 

iIa 3A5F-C3 
9.1326F-03 
2.APC7F-C2 
A.7902E-02 
7.71090-02 
l.llCAE-Cl 
1.A8A7F-01 
1.88A5F-01 
2.3023F-C1 
2.7328F-C1 
3.1721F-01 
3.617AE-01 
A.067CF-C1 
A.5193E-01 
A.9736F-01 
5.A292F-01 
5.8P57F-C1 
6.3A28F-01 
6.80C3e-Cl 
7.2P81F-C1 
7.7161E-C1 
8.17A2F-C1 
8.629AE-C1 
9.C691E-C1 
9.A755F-01 
9.P33CE-C1 
1.013CF CO 
1.0361E CO 
1.0523F CC 
1.0621F CO 
1.CA61F CO 
1.06f2E CO 
1.06C3F CC 
1.0526E CO 
1.0A31F CC 
1.032PP CO 
1.022AE CC 
1.C126F 00 
1.0CAC= 00 
9.9673F-C1 
9.9112F-01 
9.P717F-C1 
9.0A7 <;e-C 1 
9.8382E-C1 
9.8AC6F-C1 
9.8526F-01 
9.871AF-01 
9.99A7E-01 
9.92CCF-C1 
9.9A5AF-C1 
9.9693E-01 
9.99C5F-C1 
l.OOCPF CO 
1.0022F 00 


OtT 

0.0 

1.1737F-C3 

7.A721E-C3 

2.C297E-C2 

3.9193E-C2 

6.3O09E-O2 

9.084PE-02 

1.21A8F-01 

1.5A19E-01 

I. 0e37F-Ol 
2.2359E-C1 
2.5953=-Cl 
2.9597E-01 
3.3275E-C! 
3.6976E-C1 
A.0693F-01 
A.AA21F-C1 
A.3156F-C1 
S.1996F-01 
5.'5639F-C1 
5.9385F-C1 
6.3132E-C1 
6.688CE-01 
7.062eF-Cl 
7.A377E-C1 
7.8127E-01 
0.1877F-C1 
0.56CCE-C1 
8.92C5E-C1 
9.2573F-C1 
9.5599F-C1 
9.82C7F-01 
1.0035F OC 

II. 02C3F 00 
!i.C32AE CO 

1.0AC3F CC 
1.04.AAE CO 
1.CA53E CC 
1.0A36E CO 
I.OACIE CO 
1.C351E CO 
1.C29AF 00 
1.C23AF 00 
1.C17AE CO 
1.C118F 00 
1.0067E CC 
1.0C2AE CO 
9.g8eiF-ci 
9.96C7E-C1 
9.9412E-01 
9.9289E-C1 
9.9229E-01 
9.9222F-01 
9.9257F-C1 
9.9323F-C1 


Table 5.2. Computer output data for nonlinear system 
outputs and sensitivity. 
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TImf 

5.5000E 00 
5.60CCE 00 
5.7000E 00 
5.80CCE 00 

5. qOOOE 00 
6.0000E 00 
6.1000E 00 
6.200CE 00 
6.3000E 00 

6. AOOOE 00 
6.5000E 00 
6.6000E 00 
6.7000E 00 

6. e000E 00 
6.9000E 00 
7.0000E 00 
7.1000E 00 
7.20C0E 00 
7.3000E 00 
7.4000E 00 
7.5000E 00 
7.6000E 00 
7.7000E 00 

7. a000E 00 
7.90C0E 00 
8.0000E 00 
8.10C0E 00 
8.20C0E 00 
8.3000E 00 
8.4000E 00 
8.5000E 00 
8.60C0E 00 
8.7000E 00 
8.8000E 00 
8.9000E 00 
9.0000E 00 
9.10008 00 
9.2000E 00 
9.3000E OC 
Q.AOCOE 00 
<J.5000E 00 
9.60C0E 00 
9.7000E 00 
9.8000P 00 
9.9000E 00 
l.OOOOE 01 
I.OIOOE 0! 
1.0200E 01 
1.0300E 01 
1.0400F 01 
1.0500E 01 
1.0600P 01 
1.0700E 01 
1.0800E 01 
1.0900E 01 


Table 5.2. 


5.2703E-C2 
4.98778-02 
4.48898-02 
3.8392E-02 
3.10018-02 
2.3272E-02 
1.5682E-02 
8.6141E-C3 
2.3564E-03 
2.8982E-C3 
7.0448E-03 
1.0060E-02 
1.1984P-02 
1.2^08E-02 
1.29£2E-02 
1.2303E-C2 
1.1094E-C2 
9.5025E-03 
7.6836E-03 
5.7786E-03 
3.9067E-03 
2.1646E-03 
6.2460E-04 
6.6608E-04 
1.68298-03 
2.4202E-03 
2.88898-03 
3.1124E-03 
3.1238E-03 
2.9615E-03 
2.6658E-03 
2.2785E-C3 
1.8387e-C3 
1.3797E-03 
9.3007E-04 
5.1326E-04 
1.4567E-04 
1.6022E-04 
4.00598-04 
5.7411E-04 
6.8283E-04 
7.3338E-CA 
7.3433E-04 
6.9433E-04 
6.2269E-04 
5.2989E-04 
4.2456F-C4 
3.1644E-C4 
2.1058E-04 
1.1283E-C4 
2.7299E-05 
4.33338-05 
9.8586E-C5 
1.38168-04 
1.6284E-C4 


K=50 

OLT 

9.9813F-01 
<?.Qq44p_oi 
1.00C5E CO 
1.00138 00 
1.CC19F 00 
1.00238 CO 
1.0024F 00 
1.0024F 00 
1.0023F CO 
1.002CE OOi 
1.CC17E 00 
1.00138 00 
I.OCICF 00 
1.00C6F 00 
1.00038 00 
I.OOOIF 00 
9.9985F-C1 
9.99698-01 
9.9958F-C1 
9.9951E-C1 
9.9949E-01 
9.995CF-C1 
9.9953F-C1 
9.9959F-01; 
9.9966F-01 
9.9973F-01 
9.99R1F-C1 
9.9Ca8F-Cl 
9.9994F-C1 

9.99998-01 
l.OOOOE 00 
I.OOOIE 00 
I.OCCIF 00 
I.OOOIE 00 
I.OOOIF 00 
I.COCIF CO 
I.OOOIE 00 
I.OOOIE CO 
I.OOOIF CO 
I.OOOIE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE 00 
l.OOOCF 00 
9.9999F-C1 

9.99988- 01 
9.9998F-C1 

9.99989- 01 
9.9C98F-01 

9.99988-01 
9.999PE-01 

9.99988-01 

9.99988-01 

9.99998-01 
9.9999F-C1 


K=55 

OLT 

1.0031F 00 

1.0037F 00 

1.C04CF 00 
1.0C39F CO 
1.0036F CO 
1.C031E 00 
1.C026E CO 
I.OCICF 00 

1.0013E 00 
1.00078 OC 
1.00028 00 
9.99808-01 
9.9946F-01 
9.99238-01 
9.99098-01 
9.9903F-01 
9.9904F-C1 
9.9912E-01 
9.9923F-01 
9.9937P-CI 
9.9952F-C1 
9.99678-01 
9.9982F-C1 

9.99948-01 
l.OOOOE 00 
I.OOOIF 00 
1.0CC2F CO 
1.00028 CO 
1.00028 00 
1.CCC2E 00 
1.0CC2E CO 
1.C002F 00 
1.00C2E CO 
I.OOOIF 00 
I.OOOIF 00 
1.C0CC8 CO 
l.OOOOE 00 
9.9999F-01 
9.99978-01 
9.99958-01 

9.99948-01 

9.99948-01 
9.9994F-01 

9.99948-01 
9.99958-01 
9.9996E-C1 
9.9997E-C1 
9.9998F-01 

9.99998-01 
9.99C9E-C1 
l.OOOOE 00 
l.OOOOE 00 
l.OOOOE CO 
l.OOOOE CO 
l.OOOOE CO 


K=45 

OLT 

9.9411E-C1 
9.9511E-01 
9.9616F-C1 
9.9718E-C1 
9.9813E-C1 
9.9899E-01 
9.99718-01 
1.00C3E 00 
1.00C7E 00 
l.OOllF CO 
1.0012F 00 
1.0013F 00 
1.0013E 00 
1.CC13F 00 
l.COllE 00 
I.OOICE OC 
1.00C8E 00 
1.00C6F CO 
1.00C5E 00 
1.00C3E CO 
I.OOOIE 00 
l.OOOOE 00 
9.9993E-01 
9.9986E-C1 
9.99P1E-01 
9.9978F-01 
9.9977F-C1 
9.9977E-01 
9.9970E-O1 
9.998CE-C1 
9.9983E-C1 
9.9986E-C1 
9.9989E-C1 
9.9992F-01 
9.9995E-01 
9.9997F-01 
9.9999F-01 
l.OOOOE 00 
l.OOOOE CO 
l.OOOCF CO 
1.0000= CO 
l.OOOCF 00 
l.OOOOE CO 
l.OOOOE CO 
l.OOOCF 00 
l.OOOOE 00 
l.OOOCF 00 
l.CCCCE 00 
l.OOOCF 00 
l.OOOCF CO 
l.OOOCF CO 
l.OOOOE 00 
9.9999E-C1 
9.9999F-01 
9.9999E-01 


(Continued) 
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VI. CONCLUSIONS 


The technique of generating sensitivity functions for 
a basic linear feedback control system using the concept 
of a sensitivity model is developed. This model‘is a replica 
of the original system and is cascaded with it. All sensi¬ 
tivity functions are generated simultaneously with the time 
response by use of the IBM/360 Digital Simulation Program. 

Sensitivity of a differential equation with respect to 
its coefficients is investigated. Computer graphs show the 
differences in the influence of the coefficients on the 
solution of the differential equation. The coefficients of 
the lower-order terms have the greatest influence but the 
effect is obtained with the greatest time delay. 

A theoretical sensitivity model is developed for a non¬ 
linear control system. Instead of a replica of the original 
control system as required for the linear model, the non¬ 
linear element is replaced by a component with the instanta¬ 
neous gain equal the slope of the nonlinear function. An 
example illustrates the procedure for a nonlinearity due to 
a saturation element. 

Sensitivity analysis indicates that proper interpretation 
of the sensitivity function results in accurate estimates for 
changes in the system response due to parameter variations. 
Tolerance limits of parameter variations are determined and 
performance improved by application of sensitivity coefficients. 
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APPENDIX A 


DSL/360 DIGITAL SI^'Ul.ATICN PROGRAM FOR BASIC SYSTEM 


TITLE BASIC SEPVC SENSITIVITY COE'=FIC lENTS 
INTEG RKS 

STORAG IC(3ItNUM(l),DEN(4) 

TABLE IC( l-3)=4*0tMJM( l)=1.0,DEN(l-4) =1. ,10.,2^ . ,0. 
CONTRL FINTIH=IO.♦DELT=.l 
PARAM IN = 1. ,F = 1.,K=50. 

♦ SOLUTION C«= BASIC SYSTEM FCR TIME RESPONSE 

ERRORl=IK-F*CUT 

OUT=TRNFR(0»3,IC,NUM,DEN»K*ERR0R1) 

♦ SOLUTION OF MODEL FOR SENSITIVITY COEFFICIENTS 

♦ SENSK=SENSITIVITY TO FORWARD GAIN K 

♦ SENSF=SENSITIVITY TO FEEDBACK GAIN F 

SFNSK=CLT+SFNSF 

X=TRNFR(0»3,ICtNUM,DEN,K*S=NSK) 

SENSF=X*(-F) 

PRINT .1,0UT,SENSK,SFNSF 

PREPAR .1tOUTtSENSKtSENSF 

GRAPH SAME,TIME,CUTtSENSK,SENSF 

PRPLQT ONLY 

END 

STOP 
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APPENDIX P 


DSL/’dC DIGITAL SI^ULATICN FRCGRAM FOR MULTILOOP SYSTEM 


TITLE MULTILCCP SERVO SENSITIVITY COEFFICIENTS 
INTEG RKS 

CONTRL FINTIM=IC.0»DELT=.01 

PARAV Kl=40.tK2=l.♦K3 = l.,FGAIN=I.♦IN=1.♦TAU 1 = 4 . , ... 
TAU4=1. 

* SOLLTICN CF MLLTILCCP SYSTFV FCR TIME RESPONSE 

ERROR 1 = 1N-fgaIN^CUT 
0=EXCIT*(-KA)*TAUA 
F6Dl=n 

FED2=RFALPL(0.tTAU4,D) 

FEDP=FEC1-FEC2 
EPR0R2=FRRCRI+FE0B 
EXCIT=REALPL(0.tTAUl♦Kl*FRRnR2) 

GEN=REALPL(0.♦2. ,K2*5XCIT) 

0UT=PFALPL(0.,1.,K3*GFN) 

* SOLUTION CF VCTEL FOR SENSITIVITY COEFFICIENTS 

SFNSK3=CLT-fgAIN*Z 
nx=X*(-KA ) *TAIIA 
FEDXl=CX 

FEDX?=PFALPL(0. .TAUA^DX) 
cpnSKA=FEDX1-FE0X2 
SENSTA=REALPL(0.♦TAUA,S=NSK4) 

SFNSK1=SENSK3+SENSK4 

P=-SFNCK1 

0=RFALPL(0.tTAUl»P) 

SFNSTl=P-Q 

X=RFALPL< 0. ,4. ,KI*SENSK1) 

Y=REALPL(0.^2.♦K2*X) 

Z=RFALFL(0.,1.tKBAY) 

PRINT .1tOLT ,SENSK3 tSENSKl ,SENST1tSENSKA,SENST4 

PREPAR .1,0UT»SFNSK3iSFNSKl .SENS-^l .SENSKA.SENSTA 

GRAPH SAME,TIME,CUT,SENSK3,SENSK1,SENSEI»SENSK4»SFNST4 

GRAPH SAmE»TIME tSENSKltSENSTl 

PRPLOT ONLY 

END 

STOP 
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APPFNDIX C 


DSL/36C digital SI^'ULATICN PROGRAM' FOR SIXTH ORDER 
DIFFERENTIAL FQLATIOM 


TITLE SIXTH ORDER DIFFERENTIAL EQUATION SENSITIVITY MODEL 
IN^EG RKS 

CONTPL FINTIM=50.,D6LT=.05 

PARAM A = 1 . ♦B = 3.A ,C=11.4,0=9.4»E=7.4,F = .S75,G=.45 2 

* SOLUTION CF DE FCR TIME RESPONSE 

IN=STFP(C) 

XD6 = (--B*X05I+(-C*X04) + (-D*XD3) + (-E*XD2 ) + (-F*XDl ) .. . 

+(-G*X»+IN 

XD5 = IMGRL(0.tXD6) 

XD4=INTGRL(0.,XD5) 

X03=INTGRL(0.,X04) 

XD2=INTGPL(0.,X03) 

XD1=INTGPL(0.,XD2) 

X=INTGRL(0.,XDl) 

* SOLUTION OF MODEL FOR SENSITIVITY COEFFICIENTS 

SFNSG=-G*SX 

SENSF=-F*SXD1 

SENSF=-F*SXD2 

SENSD=-n*SXD3 

SENSC=-C*SXD4 

SENSB=-R*SXD5 

SXD6=SENSB+SENSC+SENS0+SFNSF+SeNSF+SENS6+X 
SXD‘5 = INTGRL( 0. ,SX06) 

SXD^=INTGRL(C.»SX05) 

SXD3=INTGPL(C.♦SXD4) 

SXD2 = IN'TGRL(0. tSXD3) 

SX01=INTGRL(C.»SXD2) 

SX=INTGRL(0.»SXDl) 

PRINT .1tXtSENSBtSENSCtSENSDtSENSE,SENSF,SENSG 

PREPAR .1,X,SFNS8,SENSC,SENSD,SENSE,SENSFtSENSG 

GRAPH SAME,TIME,X,SENSB,SENSC,SENSD,SENSE,SENSE, SENSG 

PRPLOT ONLY 

END 

STOP 
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APPENDIX D 


DSL/360 DIGITAL SIf'ULATICN PROGRAM FOR NONLINEAR STSTFM 


TITLF NONLINEAR SYSTEM SENSITIVITY 

INTEG RKS 
STORAG IC(3),8(AJ 

TABLE IC( 1-3) = 3>*Q,B( 1-4) =1. ,10. ,24. »0. 

CQNTRL FINTIM=12.,DELT=.l 
PARAM K = 50.,Pl=-.2 ,P? = .2 

* SOLUTION CF NCM.INFAP SYSTEM FOR TIME RESPONSE 

IN=STFP(C) 

ERROR=TN-CUT 
Y = LIMIT(P1 ,P2,ERROR) 

OUT=DECN(3,10,8,K+Y) 

* SOLUTION OF mqceL FOR SENSITIVITY COEFFICIENTS 

* SENSK=SENSITIVITY TO FORWARD GAIN K 

NOSORT 

IFIERRCR.GT.P2.OR.ERROR.LT.PDGC TO 1 

M=OUT-N 

GO TO 2 

1 M = C 

2 CONTINUE 
N=DEQN(3,IC,B,K*M) 

SFNSK=CLT-N 

PRINT .1,FRRCR,CLT,SENSK 

PREPAR .l,0L•r,Sc^5K 

GRAPH SAME,TIME,CUT,SENSK 

PRPLOT ONLY 

END 

STOP 
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